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ABSTRACT 


Predicting  the  liquefaction  resistance  of  soil  is  an  important  step  in  the  engineering 
design  of  new  and  the  retrofit  of  existing  structures  in  earthquake-prone  regions.  The 
procedure  currently  used  in  the  U.S.  and  throughout  much  of  the  world  to  predict  liquefaction 
resistance  is  termed  the  simplified  procedure.  This  simplified  procedure  was  originally 
developed  by  H.  B.  Seed  and  I.  M.  Idriss  in  the  late  1960s  using  blow  count  from  the  Standard 
Penetration  Test.  Small-strain,  shear  wave  velocity  measurements  provide  a promising 
alternative  and/or  supplement  to  the  penetration-based  approach.  This  report  presents  draft 
guidelines  for  evaluating  liquefaction  resistance  using  shear  wave  velocity  measurements. 
These  draft  guidelines  were  written  in  cooperation  with  industry,  researchers  and  practitioners, 
and  evolved  from  workshops  in  1996  and  1998.  The  guidelines  outline  the  development  of  a 
recommended  procedure,  which  follows  the  general  format  of  the  penetration-based  simplified 
procedure.  The  proposed  procedure  has  been  validated  through  case  history  data  from  more 
than  20  earthquakes  and  70  measurement  sites  in  soils  ranging  from  clean  fine  sand  to  sandy 
gravel  with  cobbles  to  profiles  including  silty  clay  layers.  Liquefaction  resistance  curves  were 
established  by  applying  a modified  relationship  between  shear  wave  velocity  and  cyclic  stress 
ratio  for  constant  average  cyclic  shear  strain  suggested  by  R.  Dobry.  These  curves  correctly 
predict  moderate  to  high  liquefaction  potential  for  over  95  % of  the  liquefaction  case  histories, 
and  are  shown  to  be  consistent  with  the  penetration-based  curves.  To  further  validate  the 
procedure,  additional  case  histories  are  needed  with  all  soil  types  that  have  and  have  not 
liquefied,  particularly  from  deeper  deposits  (depth  > 8 m)  and  from  denser  soils  (shear  wave 
velocity  > 200  m/s)  shaken  by  stronger  ground  motions  (peak  ground  acceleration  > 0.4  g). 
The  guidelines  serve  as  a resource  document  for  a final  recommended  practice,  and  for 
practitioners  and  researchers  involved  in  evaluating  soil  liquefaction  resistance. 


KEYWORDS:  building  technology;  earthquakes;  in  situ  measurements;  seismic  testing;  shear 
wave  velocity;  soil  liquefaction 
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CHAPTER  1 


INTRODUCTION 


1.1  BACKGROUND 

A major  cause  of  damage  from  earthquakes  is  liquefaction-induced  ground  failure.  For 
example,  direct  property  loss  caused  by  liquefaction  during  the  1 989  Loma  Prieta,  California 
earthquake  (moment  magnitude,  = 7.0)  was  over  $100  million  (Holzer,  1998).  Large 
indirect  property  loss  by  fire  almost  occurred  in  1989  when  liquefaction-induced  ground 
deformation  ruptured  water  mains  that  served  the  Marina  District  of  San  Francisco. 
Fortunately,  the  fire  in  the  Marina  District  at  Divisadero  and  Beach  Streets  was  contained  to  the 
three-story  apartment  building  where  it  ignited.  It  was  also  fortunate  that  the  1989  earthquake 
did  not  occur  closer  to  the  San  Francisco  Bay  area.  The  cities  of  Kobe  and  Osaka,  Japan  were 
not  so  fortunate.  The  1995  Hyogoken-Nanbu  earthquake  (M^  = 6.9)  directly  struck  this 
metropolitan  area,  causing  over  $100  billion  in  property  damage  (Kimura,  1996).  A significant 
portion  of  the  damage  in  Kobe  can  be  attributed  to  liquefaction-induced  ground  deformation. 
Predicting  soil  liquefaction  resistance  is  an  important  step  in  the  engineering  design  of  new  and 
the  retrofit  of  existing  structures  in  seismic-prone  regions. 

The  procedure  for  predicting  the  liquefaction  resistance  of  soils  currently  used  in  the 
United  States  and  throughout  much  of  the  world  is  termed  the  simplified  procedure.  This 
simplified  procedure  was  originally  developed  by  Seed  and  Idriss  (1971)  using  blow  count 
from  the  Standard  Penetration  Test  (SPT)  correlated  with  a parameter  representing  the  seismic 
loading  on  the  soil,  called  cyclic  stress  ratio.  Since  1971,  the  procedure  has  been  revised  and 
updated  (Seed,  1979;  Seed  and  Idriss,  1982;  Seed  et  al.,  1983;  Seed  et  al.,  1985).  Correlations 
based  on  the  Cone  Penetration  Test  (CPT),  the  Becker  Penetration  Test  (BPT),  and  shear  wave 
velocity  measurements  have  also  been  developed  by  various  investigators.  General  reviews  of 
the  simplified  procedure  are  contained  in  a report  by  the  National  Research  Council  (1985)  and 
a workshop  report  edited  by  Youd  and  Idriss  (1997). 

Small-strain  shear  wave  velocity,  V^,  measurements  provide  a promising  alternative 
and/or  supplement  to  the  penetration-based  approach.  The  use  of  as  an  index  of  liquefaction 
resistance  is  solidly  based  since  both  and  liquefaction  resistance  are  influenced  by  many  of 
the  same  factors  (e.g.,  void  ratio,  state  of  stress,  stress  history,  and  geologic  age). 
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The  in  situ  can  be  measured  by  several  seismic  tests  including  crosshole,  downhole, 
seismic  cone  penetrometer  (SCPT),  suspension  logger,  and  Spectral-Analysis-of-Surface- 
Waves  (SASW).  A review  of  these  test  methods  is  given  in  Woods  (1994).  ASTM  D-4428M- 
9 1 provides  a standard  test  method  for  crosshole  seismic  testing.  Standard  test  methods  do  not 
exist  for  the  other  seismic  tests. 

Some  advantages  of  using  are  (Dobry  et  al.,  1981;  Seed  et  al.,  1983;  Stokoe  et  al., 
1988a;  Tokimatsu  and  Uchida,  1990):  (1)  Measurements  are  possible  in  soils  that  are  hard  to 
sample,  such  as  gravelly  soils  where  penetration  tests  may  be  unreliable.  (2)  Measurements  can 
be  performed  on  small  laboratory  specimens,  allowing  direct  comparisons  between  laboratory 
and  field  behavior.  (3)  is  a basic  mechanical  property  of  soil  materials,  directly  related  to 
small-strain  shear  modulus,  G^,  by: 

G^  = pV/  (1.1) 


where 


p = the  mass  density  of  soil. 

(4)  G^,  or  is  in  turn  a required  property  in  analytical  procedures  for  estimating  dynamic 
shearing  strain  in  soil  in  earthquake  site  response  and  soil-structure  interaction  analyses.  (5) 
can  be  measured  by  the  SASW  test  method  at  sites  where  borings  may  not  be  permitted,  such 
as  capped  landfills,  and  at  sites  that  extend  for  great  distances  where  rapid  evaluation  is 
required,  such  as  lifelines  and  large  building  complexes. 

Three  concerns  when  using  to  evaluate  liquefaction  resistance  are:  (1) 
Measurements  are  made  at  small  strains,  whereas  pore-water  pressure  buildup  and  liquefaction 
are  medium-  to  high-strain  phenomenon  (Jamiolkowski  and  Lo  Presti,  1990; 
Teachavorasinskun  et  al.,  1994;  Roy  et  al.,  1996).  This  concern  can  be  significant  for  cemented 
soils,  since  small-strain  measurements  are  highly  sensitive  to  weak  interparticle  bonding  which 
is  eliminated  at  medium  and  high  strains.  It  can  also  be  significant  in  silty  soils  above  the  water 
table  where  negative  pore  water  pressures  can  increase  V^.  (2)  No  samples  are  obtained  for 
classification  of  soils  and  identification  of  non-liquefiable  soft  clayey  soils.  According  to  the 
so-called  Chinese  criteria,  non-liquefiable  clayey  soils  have  clay  contents  (particles  smaller 
than  5 pm)  greater  than  15  %,  liquid  limits  greater  than  35  %,  or  moisture  contents  less  than 
90  % of  the  liquid  limit  (Seed  and  Idriss,  1982).  (3)  Thin,  low  strata  may  not  be  detected  if 
the  measurement  interval  is  too  large  (USER,  1989;  Boulanger  et  al.,  1997). 
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In  general,  borings  should  always  be  a part  of  the  field  investigation.  Surface 
geophysical  measurements  and  cone  soundings  are  often  conducted  first  to  help  select  the  best 
locations  for  borehole  sampling  and  testing.  Surface  geophysical  tests  usually  involve  making 
measurements  at  several  different  locations,  and  provide  general,  or  average,  stratigraphy  for 
sediments  beneath  the  area  tested.  The  ability  of  surface  geophysical  methods  to  resolve  a 
layer  at  depth  depends  on  the  thickness,  depth,  and  continuity  of  that  layer,  as  well  as  the  test 
and  interpretation  procedures  employed.  Cone  soundings  provide  detailed  stratigraphy  at  each 
test  location  for  sediments  that  can  be  penetrated.  The  preferred  practice  when  using 
measurements  to  evaluate  liquefaction  resistance  is  to  drill  sufficient  boreholes  and  conduct 
sufficient  tests  to  detect  and  delineate  thin  liquefiable  strata,  to  identify  non-liquefiable  clay- 
rich  soils,  to  identify  silty  soils  above  the  ground  water  table  that  might  have  lower  values  of 
should  the  water  table  rise,  and  to  detect  liquefiable  weakly  cemented  soils. 


1.2  PURPOSE 

This  report  presents  draft  guidelines  for  evaluating  liquefaction  resistance  through  shear 
wave  velocity  measurements.  The  draft  guidelines  incorporate  suggestions  from  two 
workshops.  The  first  workshop  was  held  on  January  4-5,  1996  in  Salt  Lake  City,  Utah,  and  was 
sponsored  by  the  National  Center  for  Earthquake  Engineering  Research  (NCEER).  The  second 
workshop  was  held  on  August  14-15,  1998  also  in  Salt  Lake  City,  and  was  sponsored  by  the 
Multidisiplinary  Center  for  Earthquake  Engineering  Research  (MCEER,  formally  NCEER)  and 
the  National  Science  Foundation  (NSF).  The  guidelines  outline  the  development  of  a 
recommended  procedure  based  on  the  suggestions  given  at  these  two  workshops,  herein  called 
the  1996  NCEER  workshop  and  the  1998  MCEER  workshop.  The  guidelines  provide  guidance 
on  selecting  site  variables  and  correction  factors  that  are  consistent  with  the  shear-wave-based 
procedure. 


1.3  REPORT  OVERVIEW 

Following  this  introduction.  Chapter  2 outlines  the  development  of  several  proposed 
relationships  between  liquefaction  resistance  and  V^.  Chapter  3 presents  case  history  data  and 
describes  their  general  characteristics.  Chapter  4 establishes  the  recommended  liquefaction 
resistance  evaluation  curves  from  the  case  history  data.  Chapter  5 shows  how  the 
recommended  evaluation  curve  is  applied,  as  demonstrated  by  two  case  studies.  And  Chapter  6 
summarizes  the  recommended  procedure,  as  well  as  identifies  issues  that  remain  to  be  resolved. 
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To  assist  the  reader,  Appendix  A provides  a list  of  Symbols  and  Notation,  and 
Appendix  B provides  a Glossary  of  Terms.  Appendix  C presents  a summary  of  case  history 
data  used  to  develop  the  recommended  curves. 
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CHAPTER  2 


LIQUEFACTION  RESISTANCE  AND  SHEAR  WAVE  VELOCITY 


During  the  past  two  decades,  several  procedures  for  predicting  liquefaction  resistance 
based  on  have  been  proposed.  These  procedures  were  developed  from  laboratory  studies 
(Dobry  et  al.,  1981;  Dobry  et  al.,  1982;  de  Alba  et  al.,  1984;  Hynes,  1988;  Tokimatsu  and 
Uchida,  1990;  Tokimatsu  et  al.,  1991a;  Rashidian,  1995),  analytical  studies  (Bierschwale  and 
Stokoe,  1984;  Stokoe  et  al.,  1988c;  Andrus,  1994),  penetration- correlations  (Seed  et  al., 
1983;  Lodge,  1994;  Kayabali,  1996;  Rollins  et  al.,  1998b),  or  field  performance  data  and  in  situ 
measurements  (Robertson  et  al.,  1992;  Kayen  et  al.,  1992;  Andrus  and  Stokoe,  1997). 
Several  of  these  procedures  follow  the  general  format  of  the  simplified  procedure,  where  is 
corrected  to  a reference  overburden  stress  and  correlated  with  the  cyclic  stress,  or  resistance, 
ratio. 


2.1  CYCLIC  STRESS  RATIO 


The  cyclic  stress  ratio,  CSR,  at  a particular  depth  in  a level  soil  deposit  can  be  expressed 
as  (Seed  and  Idriss,  1971): 


r 


CSR  = — = 0.65 

cr'v 


CLr 


V S J \G  V J 


(2.1) 


where 

= the  average  equivalent  uniform  shear  stress  generated  by  the  earthquake 
assumed  to  be  0.65  of  the  maximum  induced  stress, 

^max  - horizontal  ground  surface  acceleration, 

(j'v  = the  initial  effective  vertical  (overburden)  stress  at  the  depth  in  question, 

CTv  = the  total  overburden  stress  at  the  same  depth, 

g = the  acceleration  of  gravity,  and 

= a shear  stress  reduction  coefficient  to  adjust  for  flexibility  of  the  soil  profile. 
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Equation  2.1  is  based  on  Newton’s  second  law  where  force  is  equal  to  mass  times  acceleration. 
The  coefficient  is  added  because  the  soil  column  behaves  as  a deformable  body  rather  than  a 
rigid  body. 

2.1.1  Peak  Horizontal  Ground  Surface  Acceleration 

Peak  horizontal  ground  surface  acceleration  is  a characteristic  of  the  ground  shaking 
intensity,  and  is  defined  as  the  peak  value  in  a horizontal  ground  acceleration  record  that  would 
occur  at  the  site  without  the  influence  of  excess  pore-water  pressures  or  liquefaction  that  might 
develop  (Youd  et  al.,  1997).  Peak  accelerations  are  commonly  estimated  using  empirical 
attenuation  relationships  of  as  a function  of  earthquake  magnitude,  distance  from  the 
energy  source,  and  local  site  conditions. 

Regional  or  national  seismic  hazard  maps  (Frankel  et  al.,  1996;  Frankel  et  al.,  1997; 
http://geohazards.cr.usgs.gov/eq/)  are  also  often  used  to  estimate  peak  accelerations.  If  peak 
acceleration  is  estimated  from  a map,  the  magnitude  and  distance  information  should  be 
obtained  from  the  deaggregated  matrices  used  to  develop  the  map.  The  value  of  selected 
will  depend  on  the  target  level  of  risk  and  compatibility  of  site  conditions.  For  site  conditions 
not  compatible  with  available  probabilistic  maps  or  attenuation  relationships,  the  value  of 
may  be  corrected  based  on  dynamic  site  response  analyses  or  site  class  coefficients  given  in  the 
latest  building  codes. 

2.1.2  Total  and  Effective  Overburden  Stresses 

Required  in  the  calculation  of  Gv  and  g\  are  densities  of  the  various  soil  layers,  as  well 
as  characteristics  of  the  ground  water.  For  non-critical  projects  involving  hard-to-sample  soils 
below  the  ground  water  table,  densities  are  often  estimated  from  typical  values  for  soils  with 
similar  grain  size  and  penetration  or  velocity  characteristics.  Fortunately,  CSR  is  not  very 
sensitive  to  density,  and  reasonable  estimates  of  density  yield  reasonable  results. 

The  values  of  g\  and  CSR  are  sensitive  to  the  ground  water  table  depth.  Other  ground 
water  characteristics  that  may  be  significant  to  liquefaction  evaluations  include  seasonal  and 
long-term  water  level  variations,  depth  of  and  pressure  in  artesian  zones,  and  whether  the  water 
table  is  perched  or  normal. 

2.1.3  Stress  Reduction  Coefficient 

2. 1.3.1  Relationship  by  Seed  and  Idriss  (1971) — Values  of  are  commonly 
estimated  from  the  chart  by  Seed  and  Idriss  (1971)  shown  in  Fig.  2.1.  This  chart  was 
determined  analytically  using  a variety  of  earthquake  motions  and  soil  conditions.  Average 
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Stress  Reduction  Coefficient, 


Fig.  2.1-  Relationship  Between  Stress  Reduction  Coefficient  and  Depth  Developed  by  Seed 
and  Idriss  (1971)  with  Approximate  Average  Value  Lines  from  Eq.  2.2.  (after 
Youd  et  al.,  1997) 
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values  given  in  the  chart  can  be  estimated  using  the  following  functions  (Liao  and  Whitman, 
1986;  Robertson  and  Wride,  1997): 


r^=  1.0- 0.00765  z 

for  z < 9.15  m 

{l.ldi) 

r^=  1.174- 0.0267  z 

for  9.15  m < z ^ 23  m 

(2.2b) 

= 0.744- 0.008  z 

for  23  m < z < 30  m 

(2.2c) 

where 


z = the  depth  below  the  ground  surface  in  meters. 

Figure  2.1  shows  the  average  values  approximated  by  Eq.  2.2. 

2.1.3.2  Revised  Relationship  Proposed  by  Idriss  (1998;  1999) — ^Figure  2.2  presents 
revised  average  values  of  proposed  by  Idriss  (1998;  1999)  for  various  earthquake 
magnitudes.  The  plotted  curves  are  averages  of  many  individual  curves  derived  analytically  by 
Golesorkhi  (1989)  under  the  supervision  of  the  late  Prof.  H.  B.  Seed.  They  are  defined  by  the 
following  relationship  (after  Idriss,  1998;  modified  for  depth  in  meters): 

ln{r^  = a{z)  + p(z)  (2.3) 


where 


a{z)  = -1.012  - 1.126  sin 


f 


Ul.7 


+ 5.133  , 


(2.4) 


p{z) 


= 0.106  + 0.1 18  sin 


+ 5.142  , 


(2.5) 


As  shown  in  Fig.  2.2,  the  curve  defined  by  Eq.  2.3  for  = 7.5  is  almost  identical  to  the 
average  of  the  range  published  by  Seed  and  Idriss  (1971). 

The  scatter  in  the  individual  curves  used  to  determined  the  average  curves  shown  in  Fig. 
2.2,  as  well  as  Fig.  2.1,  is  rather  large.  For  example,  coefficients  determined  for  a 30  m thick, 
loose  sand  deposit  and  magnitude  5.5  earthquakes  exhibit  standard  deviations  of  about  0.1  at  a 
depth  of  5 m and  0.15  at  a depth  of  10  m.  These  standard  deviation  values  would  be  larger  if 
soil  deposits  of  various  thicknesses  and  densities  are  considered.  Figure  2.2  provides  an 
estimate  of  the  effects  of  earthquake  magnitude  on  r^. 
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Stress  Reduction  Coefficient,  r^j 

0.0  0.2  0.4  0.6  0.8  1.0 


Fig.  2.2  - Relationship  Between  Average  Stress  Reduction  Coefficient  and  Depth  Proposed 
by  Idriss  (1998;  1999)  with  Average  of  Range  Determined  by  Seed  and  Idriss 
(1971). 
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2.2  STRESS-CORRECTED  SHEAR  WAVE  VELOCITY 


As  mentioned  in  Chapter  1,  the  in  situ  can  be  measured  by  a number  of  methods. 
The  accuracy  of  these  methods  can  be  sensitive  to  procedural  details,  soil  conditions,  and 
interpretation  techniques. 

One  important  factor  influencing  is  state  of  stress  in  soil  (Hardin  and  Drnevich, 
1972;  Seed  et  al.,  1986).  Laboratory  tests  (Roesler,  1979;  Yu  and  Richart,  1984;  Stokoe  et  al., 
1985;  S.  H.  Lee,  1986;  N.  J.  Lee,  1993)  have  shown  that  the  velocity  of  a propagating  shear 
wave  depends  equally  on  principal  stresses  in  the  direction  of  wave  propagation,  and  the 
direction  of  particle  motion.  Thus,  measurements  made  with  wave  propagation  or  particle 
motion  in  the  vertical  direction  can  be  generally  related  by  the  following  empirical  relationship 
(Stokoe  et  al.,  1985): 


v,=A  (c7'vr((T'.r  (2.6) 

where 

A = a parameter  that  depends  on  the  soil  structure, 

(j'h  = the  initial  effective  horizontal  stress  at  the  depth  in  question,  and 

m = a stress  exponent  with  a value  of  about  0. 125. 


Following  the  traditional  procedures  for  correcting  standard  and  cone  penetration 
resistances  (Marcuson  and  Bieganousky,  1977;  Seed,  1979;  Liao  and  Whitman,  1986;  Olsen, 
1997;  Robertson  and  Wride,  1997;  Youd  et  al.,  1997;  Robertson  and  Wride,  1998),  one  can 
correct  V^to  a reference  overburden  stress  by  (Sykora,  1987b;  Robertson  et  al.,  1992): 


0.25 


where 


(2.7) 


Pa  = a reference  stress,  100  kPa  or  approximately  atmospheric  pressure,  and 

(j'v  = the  initial  effective  vertical  stress  in  kPa. 


Equation  2.7  assumes  that  cr7  = K^a'v  and  is  a constant  (~  0.5  at  sites  where  liquefaction 
has  occurred),  from  the  relationship  given  in  Eq.  2.6.  Also,  Eq.  2.7  implicitly  assumes  that 
is  measured  with  both  the  directions  of  particle  motion  and  wave  propagation  polarized  along 
principal  stress  directions  and  one  of  these  directions  is  vertical. 
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Since  the  direction  of  wave  propagation  and  the  direction  of  particle  motion  is  different 
with  respect  to  the  stress  in  the  soil  for  each  in  situ  seismic  test  method,  some  variations 
between  measured  is  expected.  These  variations  are  minimized  by  performing  the  tests  with 
at  least  a major  component  of  wave  propagation  or  particle  motion  in  the  vertical  direction.  To 
have  a major  component  of  wave  propagation  or  particle  motion  in  the  vertical  direction, 
crosshole  tests  are  conducted  with  particle  motion  in  the  vertical  direction,  and  downhole  and 
seismic  cone  tests  are  conducted  at  depths  greater  than  the  distance  between  the  source  and  the 
borehole  or  cone  sounding  such  that  wave  propagation  is  in  the  vertical  direction. 


2.3  CYCLIC  SHEAR  STRAIN 


Liquefaction  results  from  the  rearranging  of  soil  particles  and  the  tendency  for  decrease 
in  volume.  Experimental  and  theoretical  studies  show  that  decrease  in  volume  is  more  closely 
related  to  cyclic  strain  than  cyclic  stress  (Silver  and  Seed,  1971);  a threshold  cyclic  strain  exists 
below  which  neither  rearrangement  of  soil  particles  nor  decrease  in  volume  take  place 
(Dmevich  and  Richart,  1970;  Youd,  1972;  Pyke  et  al.,  1975),  and  no  pore  water  pressure 
buildup  occurs  (Dobry  et  al.,  1981;  Seed  et  al.,  1983);  and  that  there  is  a predictable  correlation 
between  cyclic  shear  strain  and  pore  pressure  buildup  of  saturated  soils  (Martin  et  al.,  1975; 
Park  and  Silver,  1975;  Finn  and  Bhatia,  1981;  Dobry  et  al.,  1982;  Hynes,  1988).  The  threshold 
cyclic  strain  is  limited  to  a narrow  range  of  variation,  ranging  from  about  0.005  % for  gravels  to 
0.01  % for  normally  consolidated  clean  sands  and  silty  sands  to  0.03  % for  overconsolidated 
clean  sands.  In  addition,  cyclic  strain-controlled  test  results  are  less  affected  than  stress- 
controlled  tests  by  factors  such  as  density,  confining  stress,  anisotropic  confining  stress,  fabric 
and  prestaining  (Martin  et  al.,  1975;  Dobry  and  Ladd,  1980;  Dobry  et  al.,  1982;  Hynes,  1988). 
It  should  also  be  noted  that  the  steady  state  approach  to  liquefaction  evaluation  by  Poulos  et  al. 
(1985)  is  based  on  a triggering  strain  level.  These  findings  confirm  the  fact  that  cyclic  strain  is 
more  fundamentally  related  to  pore  pressure  buildup  than  cyclic  stress,  and  are  strong 
arguments  in  favor  of  a cyclic  strain  approach  to  liquefaction  evaluation. 


Cyclic  shear  strain  and  cyclic  shear  stress  can  be  related  by  the  following  equation: 


(2.8) 


where 

= the  average  peak  cyclic  shear  strain  during  a cyclic  stress-controlled  test  of 
uniform  cyclic  shear  stress  which  results  in  triggering  of  liquefaction,  and 
(G)^  = the  secant  shear  modulus  at  during  the  same  cyclic  test. 
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In  the  cyclic  strain  approach  proposed  by  Dobry  et  al.  (1982),  the  average  cyclic  shear 
strain  caused  by  an  earthquake  is  estimated  from: 


= 0.65 


dmax 

g 


CvrdG  max 

pVs\G) 


(2.9) 


Equation  2.9  is  obtained  by  combining  Eqs.  1.1,  2.1  and  2.8.  The  variation  of  shear  modulus 
with  strain  is  commonly  expressed  in  terms  of  (G)^  /G^,  called  the  modulus  reduction  factor. 
The  modulus  reduction  factor  can  be  estimated  from  an  experimentally  determined  correlation. 
Neither  pore  pressure  buildup  nor  liquefaction  will  occur  when  is  less  than  the  threshold 
strain.  When  y^^  is  greater  than  the  threshold  strain,  then  pore  pressure  buildup  can  occur.  The 
amount  of  pore  pressure  buildup  can  also  be  estimated  from  an  experimentally  determined 
correlation. 


R.  Dobry  (personal  communication  to  R.  D.  Andrus,  1996)  also  derived  a relationship 
between  V^j  and  CSR  for  constant  average  cyclic  shear  strain  using  Eqs.  1.1  and  2.8. 
Combining  Eqs.  1.1  and  2.8,  and  dividing  both  sides  by  a'v  leads  to: 


(2.10) 


For  an  overburden  stress  of  100  kPa,  and  curves  of  constant  average  cyclic  strain  can 

be  expressed  by: 


CSR  = 


^ =f(yj  Vsr 

CTv 


(2.11) 


where 


AyJ  = 


yj  ^ 


MGX 


PaJ  G 


(2.12) 


Equation  2. 1 1 provides  an  analytical  basis  for  extending  liquefaction  resistance  curves  to  zero 
at  V^j  = 0,  and  provides  a means  for  establishing  curves  at  low  values  of  (say  V^j  < 125  m/s). 
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2.4  MAGNITUDE  SCALING  FACTOR 


In  developing  the  simplified  procedure.  Seed  and  Idriss  (1982)  collected  SPT  blow 
count  measurements  from  several  sites  where  surface  effects  of  liquefaction  were  or  were  not 
observed  during  earthquakes  with  magnitudes  of  about  7.5.  They  plotted  cyclic  stress  ratios 
and  corrected  blow  counts  for  the  clean  sand  (silt  and  clay  content  < 5 %)  sites,  and  drew  a 
curve  to  bound  the  liquefaction  data  points.  For  earthquakes  with  magnitude  other  than  7.5, 
Seed  and  Idriss  proposed  magnitude  scaling  factors  to  adjust  the  curve  bounding  the 
liquefaction  data  points  for  magnitude  7.5  earthquakes. 

Table  2.1  presents  the  magnitude  scaling  factors  developed  by  Seed  and  Idriss  (1982), 
as  well  as  the  magnitude  scaling  factors  developed  by  other  investigators  in  recent  years.  These 
magnitude  scaling  factors  were  derived  from  laboratory  test  results  and  representative  cycles  of 
loading  (Seed  and  Idriss,  1982;  Idriss,  personal  communication  to  T.  L.  Youd,  1995;  Idriss, 
1998;  Idriss,  1999),  correlations  of  field  performance  data  and  blow  count  measurements 
(Ambrasey,  1988;  Youd  and  Noble,  1997),  estimates  of  seismic  energy  for  laboratory  and  field 
data  (Arango,  1996),  and  correlations  of  field  performance  data  and  in  situ  measurements 
(Andrus  and  Stokoe,  1997).  Figure  2.3  shows  a plot  of  the  various  magnitude  scaling  factors 
along  with  the  range  recommended  by  the  1996  NCEER  workshop  (Youd  et  al.,  1997). 


Table  2.1  - Magnitude  Scaling  Factors  Obtained  by  Various  Investigators,  (modified  from 
Youd  and  Noble,  1997) 


Magnitude  Scaling  Factor  (MSF) 


Moment 

Magnitude, 

(1) 

Seed 

& 

Idriss 

(1982) 

(2) 

Idriss 
(personal 
commu- 
nication 
to  T.  L. 
Youd, 
1995) 

(3) 

Idriss 

(1998) 

(4) 

Idriss 

(1999) 

(5) 

Ambraseys 

(1988) 

(6) 

Youd  & Noble 
(1997) 

Pu  % 

<20  <30  <50 

(7)  (8)  (9) 

Arango 

(1996) 

(10)  (11) 

Andrus 

& 

Stokoe 

(1997) 

(12) 

5.5 

1.43 

2.20 

1.625 

1.68 

2.86 

2.86 

3.42 

4.44 

3.00 

2.20 

2.8* 

6.0 

1.32 

1.76 

1.48 

1.48 

2.20 

1.93 

2.35 

2.92 

2.00 

1.65 

2.1 

6.5 

1.19 

1.44 

1.28 

1.30 

1.69 

1.34 

1.66 

1.99 

1.60 

1.40 

1.6 

7.0 

1.08 

1.19 

1.12 

1.14 

1.30 

1.00 

1.20 

1.39 

1.25 

1.10 

1.25 

7.5 

1.00 

1.00 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0 

8.0 

0.94 

0.84 

0.88 

0.87 

0.67 

0.73 

0.75 

0.85 

0.8* 

8.5 

0.89 

0.72 

0.79 

0.76 

0.44 

0.56 

0.65* 

*Extrapolated  from  scaling  factors  for  M^.  = 6,  6.5,  7,  and  7.5  using  MSF  = 
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Fig.  2.3  - Magnitude  Scaling  Factors  Derived  by  Various  Investigators  with  Range 
Recommended  by  the  1996  NCEER  Workshop  (after  Youd  et  al.,  1997; 
Youd  and  Noble,  1997). 
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Although  the  1996  NCEER  workshop  (Youd  et  ah,  1997)  recommended  a range  of 
magnitude  scaling  factors  for  engineering  practice,  a consensus  has  not  yet  been  reached  by  the 
workshop  participants.  At  the  August  1998  MCEER  workshop,  some  concerns  were  expressed 
about  the  upper  limit  of  the  recommended  range.  Also,  a revised  set  of  magnitude  scaling 
factors  and  stress  reduction  coefficients  (see  Section  2. 1.3.2)  were  proposed  by  I.  M.  Idriss. 
The  magnitude  scaling  factors  recommended  by  the  1996  NCEER  workshop  and  the  revised 
factors  proposed  by  Idriss  (1998;  1999)  are  discussed  below. 


2.4.1  Factors  Recommended  by  1996  NCEER  Workshop 


The  magnitude  scaling  factors  recommended  by  the  1996  NCEER  workshop  (Youd  et 
al.,  1997)  can  be  represented  by: 


MSF  = 

(M.X 

I7.5J 

where 

MSF 

= the  magnitude  scaling  factor. 

K 

= moment  magnitude,  and 

n 

= an  exponent. 

(2.13) 


The  lower  bound  for  the  range  of  magnitude  scaling  factors  recommended  by  the  1 996  NCEER 
workshop  is  defined  with  n = -2.56  (Idriss,  personal  communication  to  T.  L.  Youd,  1995)  for 
earthquakes  with  magnitude  < 7.5.  The  upper  bound  of  the  recommended  range  is  defined  with 
n = -3.3  (Andrus  and  Stokoe,  1997)  for  earthquakes  with  magnitude  < 7.5.  For  earthquakes 
with  magnitude  > 7.5,  the  recommended  factors  are  defined  with  n = -2.56.  Magnitude  scaling 
factors  defined  by  Eq.  2.13  should  be  used  with  values  given  in  Fig.  2.1. 


2.4.2  Revised  Factors  Proposed  by  Idriss  (1998;  1999) 


The  magnitude  scaling  factors  proposed  by  Idriss  (1998;  1999)  are  derived  using 
laboratory  data  from  Yoshimi  et  al.  (1984)  and  a revised  relationship  between  representative 
cycles  of  loading  and  earthquake  magnitude.  The  1998  factors  are  defined  by  the  following 
equation: 

MSF  = 37.9  (Mw)'’’^'  for  > 5.75  (2. 14a) 

M5F=  1.625  forM^<5.75  (2.14b) 
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The  1999  factors  are  defined  by  the  following  equation: 


MSF=  6.9  exp|^-^l  - 0.06 
MSF=  1.82 


forM^>5.2 


forM^  <5.2 


(2.15b) 


(2.15a) 


Figure  2.3  shows  the  magnitude  scaling  factors  defined  by  Eqs.  2.14  and  2.15.  The  difference 
between  the  1998  and  1999  magnitude  scaling  factors  proposed  by  Idriss  is  small.  Magnitude 
scaling  factors  defined  by  Eqs.  2.14  and  2.15  should  be  used  with  values  given  in  Fig.  2.2. 

2.4.3  Comparison  of  Magnitude  Scaling  Factors 

The  proposed  relationships  for  MSF  can  be  compared  directly  by  combining  them  with 
the  appropriate  stress  reduction  coefficient  into  one  factor.  This  factor  is  the  product  of  and 
the  reciprocal  of  MSF.  Figure  2.4  presents  values  of  rJMSF  for  the  range  recommended  by  the 
1996  NCEER  workshop  (Youd  et  al.,  1997)  and  those  proposed  by  Idriss  (1999).  As  shown  in 
the  figure,  there  is  not  much  difference  between  the  two  sets  of  rJMSF  values  for  magnitude  of 

7.5  and  depth  less  than  11m.  At  magnitudes  near  5.5  and  shallow  depths,  the  difference 
between  rJMSF  values  proposed  by  Idriss  (1999)  and  values  recommended  by  the  1996 
NCEER  workshop  is  as  much  as  50  %. 

The  magnitude  scaling  factors  recommended  by  the  1996  NCEER  workshop  (Youd  et 
al.,  1997)  and  the  revised  magnitude  scaling  factors  proposed  by  Idriss  (1999)  will  be 
considered  in  Chapter  4 to  establish  the  recommended  liquefaction  resistance  relationship  for 
magnitude  7.5  earthquakes. 

2.5  CYCLIC  RESISTANCE  RATIO 

The  value  of  CSR  separating  liquefaction  and  non-liquefaction  occurrences  for  a given 
is  called  the  cyclic  resistance  ratio,  CRR.  Seven  proposed  relationships  between  V^j  and 
CRR  are  compared  in  Fig.  2.5  and  briefly  discussed  below. 
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Average  Stress  Reduction  Coefficient 
Magnitude  Scaling  Factor 

0 0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 


Fig.  2.4  - Variation  of  rJMSF  with  Depth  for  Various  Magnitudes  and  Proposed  Relationships. 
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Cyclic  Resistance  Ratio,  CRR 


0.6 


0.4 


0.2 


0.0 


'Curve  adjusted  using  scaling 
factor  of  1 .19  for  magnitude 
7 earthquakes 

"Approximate  curve  for  clean 
sand  & 1 5 cycles  of  loading, 
assuming  = 0.65, 

Kq  = 0.5,  = 0.9 


Tokimatsu  & 
Uchida  (1990) 


Liquefaction 


Robertson  et  al. 
(1992) 


Mm/ =7.5 


Lodge 

(1994)* 


Andrus  & 
Stokoe 
(1997)  • 


Kayen  et  al. 
(1992)* 


No 

Liquefaction 


0 100  200  300 

Overburden  Stress-Corrected  Shear  Wave 
Velocity,  Vsh  m/s 


Fig.  2.5  - Comparison  of  Seven  Relationships  Between  Liquefaction  Resistance 
and  Overburden  Stress-Corrected  Shear  Wave  Velocity  for  Clean 
Granular  Soils. 
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2.5.1  Relationship  by  Tokimatsu  and  Uchida  (1990) 

The  “best-fit”  curve  by  Tokimatsu  and  Uchida  (1990)  shown  in  Fig.  2.5  was  determined 
from  laboratory  cyclic  triaxial  test  results  for  various  sands  with  less  than  10  % fines  (silt  and 
clay)  and  15  cycles  of  loading.  Figure  2.6  presents  the  cyclic  triaxial  test  results.  The  solid 
symbols  in  Fig.  2.6  correspond  to  specimens  obtained  by  the  in  situ  freezing  technique.  The 
open  symbols  correspond  to  specimens  reconstituted  in  the  laboratory.  Tokimatsu  and  Uchida 
defined  the  cyclic  resistance  ratio  for  cyclic  triaxial  tests,  CRR^,  as  the  ratio  of  cyclic  deviator 
stress  to  initial  effective  confining  stress,  Odl'^o'o,  when  the  double-amplitude  (or  peak-to- 
peak)  axial  strain,  DA,  reaches  5 %.  They  measured  the  elastic  shear  modulus  of  the  specimen 
at  a shear  strain  of  10'^  % just  prior  to  the  liquefaction  test.  This  small-strain  shear  modulus 
was  normalized  to  correct  for  the  influence  of  confining  pressure  and  void  ratio  by: 


G,= 


f m) 


2/3 


(2.16) 


and 


{IM-errunf 

^ € min 


where 


(2.17) 


= the  normalized  shear  modulus, 

= the  minimum  void  ratio  determined  by  standard  test  method,  and 

TTllTl  J ^ 

a'm  - the  mean  effective  confining  stress. 


Tokimatsu  and  Uchida  selected  an  exponent  of  2/3  rather  than  1/2,  as  determined  by  Hardin 
and  Dmevich  (1972),  because  it  seemed  that  a slightly  better  correlation  could  be  obtained. 
Values  of  ranged  from  0.61  to  0.91  for  the  sands  tested.  The  actual  values  of  void  ratio  in 
each  test  were  greater  than  with  values  ranging  from  about  0.65  to  about  1.4. 


By  combining  Eqs.  1.1  and  2.16,  one  obtains  the  following  relationship  for  converting 
to  mean  stress-corrected  V^: 


V =V 

* Slm  *^5 


( 1 ^ 

0.33 

GNfiemin) 

V.  CJ  m y 

[ P J 

0.5 


(2.18) 
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Stress  Ratio  to  cause  DA  = 5 % in  1 5 cycles 


1.5 
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0.5 
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• A # ▼ Intact,  in  s/ft;  freezing 
O A V □ 0 Reconstituted 


_ 

• G Niigata  Sand 

A A Meike  Sand 

♦ Ohgishima  Sand 

T V Silica  Sand 

□ Toyoura  Sand 

0 Makuhari  Sand 

Liquefaction 


"Best  fit" 
Curve 
by 

Tokimatsu 

and 

Uchida 

(1990) 


Y Lower  bound- 

(This  Report) 


No 

Liquefaction 


200  400  600  800  1000 

Normalized  Shear  Modulus,  G/v=  Gma)/W^min)  ( 


Fig.  2.6  - Relationship  Between  Liquefaction  Resistance  and  Normalized  Shear  Modulus 
for  Various  Sands  with  Less  than  10  % Fines  Determined  by  Cyclic  Triaxial 
Testing,  (modified  from  Tokimatsu  and  Uchida,  1990) 
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where 


- mean  stress-corrected  V^,  and 

(j'm  = the  mean  effective  confining  stress  in  kgf/cm^  (1  kgf/cm^  = 98  kPa). 


Tokimatsu  and  Uchida  (1990)  suggested  using  0.65  as  an  average  value  of  for  clean  sands. 
The  overburden  stress-corrected  and  can  be  related  by: 


V =V 

^ Slm 


( 1 


KOv) 


f 


0.33 


- y, 


SI 


f— 1 

Icr'vj 


0.08 


f 


0.33 


V^'lKo) 


(2.19) 


where 


= the  coefficient  of  lateral  earth  pressure  at  rest  {=a'hlG\)- 

Values  of  for  the  “best  fit”  curve  by  Tokimatsu  and  Uchida  (1990)  shown  in  Fig.  2.5  are 
determined  (This  report;  after  Tokimatsu  et  al.,  1991a)  from  Fig.  2.6  using  Eqs.  2.18  and  2.19, 
and  assuming  - 0.5,  = 0.65,  a'm  = 100  kPa,  and  soil  density  of  1.9  Mg/m^. 

For  converting  CRR^^  to  an  equivalent  field  CRR,  Tokimatsu  and  Uchida  (1990) 
suggested  the  following  expression  originally  proposed  by  Seed  (1979): 

CRR  = ^ (CRRJ  (2.20) 


where 


= a constant  to  account  for  the  effects  of  multi-directional  shaking  with  a value 
between  0.9  and  1.0. 

Values  of  CRR  for  the  “best  fit”  curve  by  Tokimatsu  and  Uchida  shown  in  Fig.  2.5  are 
determined  from  Fig.  2.6  using  Eq.  2.20  and  assuming  = 0.5  and  = 0.9. 

Since  the  other  liquefaction  resistance  relationships  shown  in  Fig.  2.5  were  drawn  to 
bound  liquefaction  case  histories,  the  more  conservative  “lower  bound”  curve  for  the  laboratory 
test  results  by  Tokimatsu  and  Uchida  (1990)  is  also  shown.  This  curve  was  drawn  (This  report) 
from  Fig.  2.6  following  the  procedure  outlined  above. 
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2.5.2  Relationship  by  Robertson  et  al.  (1992) 

The  bounding  curve  by  Robertson  et  al.  (1992)  was  developed  using  field  performance 
data  from  primarily  sites  in  the  Imperial  Valley,  California,  along  with  data  from  four  other 
sites,  as  shown  in  Fig.  2.7.  The  soil  at  these  sites  contained  as  much  as  35  % fines.  Robertson 
et  al.  corrected  using  Eq.  2.7.  The  shape  of  their  relationship  was  based  on  the  analytical 
results  of  Bierschwale  and  Stokoe  (1984).  They  reasoned  that  the  curve  should  pass  close  to 
the  Imperial  Valley  (Wildlife  site)  data  point,  since  liquefaction  did  and  did  not  occur  at  this 
site  during  the  1987  Superstition  Hills  (M^  = 6.5)  and  Elmore  Ranch  (M^  = 6.2)  earthquakes, 
respectively.  Robertson  et  al.  used  the  magnitude  scaling  factors  suggested  by  Seed  (1979), 
similar  to  factors  listed  in  Column  2 of  Table  2.1,  to  position  their  curve  for  magnitude  7.5 
earthquakes. 

2.5.3  Relationship  by  Kayen  et  al.  (1992) 

Kayen  et  al.  (1992)  studied  four  sites  that  did  and  did  not  liquefy  during  the  1989  Loma 
Prieta,  California,  earthquake  (M^  = 7.0).  The  four  sites  are:  Port  of  Richmond,  Bay  Bridge 
Toll  Plaza,  Port  of  Oakland,  and  Alameda  Bay  Farm  Island  South  Loop  Road.  The  fines 
content  for  soils  at  these  sites  ranged  from  less  than  5 % to  as  much  as  57  %.  Values  of  were 
measured  by  the  SCPT  method  and  corrected  for  overburden  stress  using  Eq.  2.7.  Figure  2.8 
presents  their  data  and  bounding  curve.  The  curve  was  adjusted  for  magnitude  7.5  earthquakes 
assuming  a MSF  of  1.19  (see  Column  3 of  Table  2.1),  as  shown  in  Fig.  2.5. 

2.5.4  Relationship  by  Lodge  (1994) 

Lodge  (1994)  considered  the  same  sites  that  Kayen  et  al.  (1992)  studied,  as  well  as  other 
sites  shaken  by  the  1989  Loma  Prieta  earthquake.  The  curve  by  Lodge  was  developed  as 
follows.  First,  cyclic  stress  ratios  for  the  entire  soil  profile  at  each  site  were  calculated. 
Second,  available  SPT  blow  counts  were  corrected  for  overburden  pressure  and  energy.  Soil 
layers  with  high  and  low  liquefaction  potential  were  identified  with  the  procedure  of  Seed  et  al. 
(1985).  Soil  layers  with  corrected  blow  count  within  3 of  the  SPT-based  curve  were  eliminated 
due  to  uncertainties  in  the  correlation.  Third,  measurements  from  SCPT  and  crosshole  tests 
were  corrected  for  overburden  stress  using  Eq.  2.7.  Fourth,  on  a “meter  by  meter”  basis,  values 
of  and  cyclic  stress  ratio  were  plotted  for  both  layer  types,  those  which  were  predicted 
liquefiable  and  those  which  were  predicted  non-liquefiable.  Fifth,  published  data  for  sites 
shaken  by  the  1983  Borah  Peak,  Idaho,  and  1964  Niigata,  Japan,  earthquakes  were  added  to  the 
plot.  Finally,  a curve  was  drawn  to  include  all  liquefiable  layers,  as  shown  in  Fig.  2.9.  Figure 
2.5  shows  the  curve  by  Lodge  adjusted  for  magnitude  7.5  earthquakes  assuming  a MSF  of  1.19 
(see  Column  3 of  Table  2.1). 


22 
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Velocity,  Vsi,  m/s 


Fig.  2.7  - Liquefaction  Resistance  Relationship  for  Magnitude  7.5  Earthquakes  and 
Case  History  Data  from  Robertson  et  al.  (1992). 
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Fig.  2.8  - Liquefaction  Resistance  Relationship  for  Magnitude  7 Earthquake  and  Case 
History  Data  from  Kayen  et  al.  (1992). 
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Fig.  2.9  - Liquefaction  Resistance  Relationship  for  Magnitude  7 Earthquakes  and  Case 
History  Data  from  Lodge  (1994). 
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2.5.5  Relationship  by  Andrus  and  Stokoe  (1997) 


The  curve  by  Andrus  and  Stokoe  (1997)  shown  in  Fig.  2.5  was  developed  for  the 
proceedings  of  the  1996  NCEER  workshop.  Several  suggestions  were  offered  at,  and  after,  the 
workshop  concerning  how  site  variables  should  be  define,  as  well  as  the  shape  of  the  boundary 
curve  separating  liquefaction  and  no  liquefaction.  Following  the  suggestions  and  using  field 
performance  data  from  20  earthquakes  and  in  situ  measurements  from  over  50  sites  in  soils 
ranging  from  clean  fine  sand  to  sandy  gravel  with  cobbles  to  profiles  including  silty  clay  layers, 
Andrus  and  Stokoe  constructed  curves  for  uncemented,  Holocene-age  soils  with  various  fines 
content,  FC.  The  values  of  were  corrected  using  Eq.  2.7.  The  curve  for  FC  < 5 % by 
Andrus  and  Stokoe  along  with  the  case  history  data  are  presented  in  Fig.  2.10 

The  shape  of  the  curve  by  Andrus  and  Stokoe  (1997)  is  based  on  a modified  relationship 
between  V^j  and  CSR  for  constant  average  cyclic  shear  strain  suggested  by  R.  Dobry  (see 
Section  2.3).  Andrus  and  Stokoe  reasoned  that  the  curve  separating  liquefiable  and  non- 
liquefiable  soils  would  become  asymptotic  to  some  limiting  upper  value  of  V^j.  This 
assumption  is  equivalent  to  the  assumption  commonly  made  in  the  SPT-  and  CPT-based 
procedures  where  liquefaction  is  considered  not  possible  above  a corrected  blow  count  of  about 
30  (Seed  et  al.,  1985;  Youd  et  al.,  1997)  and  a corrected  tip  resistance  of  about  160  (Youd  et  al., 
1997;  Robertson  and  Wride,  1998).  Upper  limits  for  blow  count  and  V^j  are  explained  by  the 
tendency  of  dense  soils  to  exhibit  dilative  behavior  at  large  strains,  causing  negative  pore  water 
pressures.  While  it  is  possible  in  a dense  soil  to  generate  pore  water  pressures  close  to  the 
confining  stress  if  large  cyclic  strains  or  many  cycles  are  applied  to  the  soil,  the  amount  of 
water  expelled  during  reconsolidation  is  dramatically  less  for  dense  soils  than  for  loose  soils. 
As  explained  by  Dobry  (1989),  in  dense  soils,  settlement  is  insignificant  and  no  sand  boils  or 
engineering  failure  take  place  because  of  the  small  amount  of  water  expelled.  This  is  important 
because  the  definition  of  liquefaction  used  to  classify  the  case  histories  here,  as  well  as  in  the 
penetration-based  simplified  procedures,  is  based  on  surface  manifestations. 

Thus,  Andrus  and  Stokoe  (1997)  modified  Eq.  2.1 1 to: 


(2.21) 


where 


V*si  = the  limiting  upper  value  of  for  liquefaction  occurrence,  and 
a,b  = curve  fitting  parameters. 
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Fig.  2.10  - Liquefaction  Resistance  Relationship  for  Magnitude  7.5  Earthquakes  and 
Uncemented  Clean  Soils  of  Holocene  Age  with  Case  History  Data  from 
Andrus  and  Stokoe  (1997). 
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The  first  term  in  Eq.  2.21  is  a form  of  Eq.  2.11,  assuming is  independent  of  initial 
effective  confining  pressure  and  of  pore  water  pressure  buildup.  The  second  term  is  a 
hyperbola  with  a small  value  at  low  values  of  and  a very  large  value  as  approaches  Vsi  • 

The  curve  by  Andrus  and  Stokoe  (1997)  shown  in  Figs.  2.5  and  2.10  is  defined  by  Eqs. 
2.13  and  2.21  with  a = 0.03,  b = 0.9,  n = -3.3,  andyji  = 220  m/s. 

2.5.6  Relationship  Proposed  in  This  Report 

Since  the  publication  of  the  1996  NCEER  workshop  proceedings  (Youd  and  Idriss,  eds., 
1997),  the  case  history  data  compiled  by  Andrus  and  Stokoe  (1997)  have  been  revised,  based 
on  new  information,  and  expanded  to  include  field  performance  data  from  26  earthquakes  and 
more  than  70  measurements  sites.  Also,  the  1998  MCEER  workshop  was  held  to  discuss 
developments  since  the  1996  workshop.  From  the  suggestions  given  at  the  second  workshop 
and  using  the  expanded  database,  the  curve  proposed  by  Andrus  and  Stokoe  (1997)  is  revised  in 
this  report.  The  case  history  data  and  the  revised  curve  for  uncemented  soils  with  FC  < 5 % are 
shown  in  Fig.  2.1 1.  Chapter  4 discusses  the  development  of  the  revised  curve. 


2.6  SUMMARY 

A simplified  procedure  for  evaluating  liquefaction  resistance  of  soils  using 
measurements  was  outlined  in  this  chapter.  Also  discussed  are  seven  proposed  relationships 
between  CRR  and  V^;.  Many  of  the  differences  among  the  seven  curves  (see  Fig.  2.5)  can  be 
explained  by  the  following  three  factors:  (1)  The  “best-fit”  curve  by  Tokimatsu  and  Uchida 
(1990)  is  more  of  a median  curve,  while  the  other  curves  bound  the  liquefaction  case  history 
data.  (2)  Portions  of  the  proposed  curves  are  based  on  limited  data,  and  the  investigator(s)  have 
assumed  different  levels  of  conservatism.  (3)  Methods  for  selecting  some  site  variables  and 
correction  factors  are  different  among  investigator(s). 
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Fig.  2.1 1 - Revised  Liquefaction  Resistance  Relationship  for  Magnitude  7.5  Earthquakes 
and  Uncemented  Clean  Soils  of  Holocene  Age  with  Case  History  Data  from 
This  Report. 
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CHAPTER  3 


CASE  HISTORY  DATA  AND  THEIR  CHARACTERISTICS 


Shear  wave  velocity  measurements  have  been  made  for  field  liquefaction  studies  at 
many  sites  during  the  past  fifteen  years.  Table  3.1  lists  over  70  sites  and  26  earthquakes  that 
have  been  investigated.  Of  the  26  earthquakes  listed,  9 occurred  in  the  United  States;  and  the 
other  15  in  Japan,  Taiwan,  and  China.  The  field  performance  information  for  these  earthquakes 
along  with  the  measurements  provides  an  important  opportunity  to  determine  the 
relationship  between  liquefaction  resistance  and  directly  from  case  histories.  A summary  of 
available  case  history  data  is  presented  in  Appendix  C.  This  chapter  describes  the  site  variables 
and  characteristics  of  the  database. 


3.1  SITE  VARIABLES  AND  DATABASE  CHARACTERISTICS 

3.1.1  Earthquake  Magnitude 

Earthquake  magnitudes  for  the  26  earthquakes  listed  in  Table  3.1  range  from  5.3  to  8.3, 
based  on  the  moment  magnitude  scale.  Moment  magnitude  is  the  scale  most  commonly  used 
for  engineering  applications,  and  is  the  preferred  scale  for  liquefaction  resistance  calculations 
(Youd  et  al.,  1997).  When  other  magnitude  scales  are  reported  by  the  investigator(s),  they  are 
converted  to  using  the  relationship  of  Heaton  et  al.  (1982)  shown  in  Fig.  3.1. 

3.1.2  Shear  Wave  Velocity  Measurement 

Shear  wave  velocity  measurements  were  made  with  139  test  arrays  at  the  more  than  70 
investigation  sites  listed  in  Table  3.1.  A test  array  is  defined  in  this  report  as  the  two  boreholes 
used  for  crosshole  measurements,  the  borehole  and  source  used  for  downhole  measurements, 
the  cone  sounding  and  source  used  for  seismic  cone  measurements,  the  borehole  used  for 
suspension  logging  measurements,  or  the  line  of  receivers  used  for  Spectral-Analysis-of- 
Surface-Waves  (SASW)  measurements.  Of  the  139  test  arrays,  39  are  crosshole,  21  downhole, 
27  seismic  cone,  15  suspension  logger,  36  SASW,  and  one  is  unknown. 
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Table  3.1  - Earthquakes  and  Sites  Used  to  Establish  Liquefaction  Resistance  Curves 


Earthquake 

Moment 

Magnitude 

Site 

Reference 

(1) 

(2) 

(3) 

(4) 

1906  San  Francisco,  Calif. 

7.7 

Coyote  Creek;  Salinas  River 

Youd  & Hoose  (1978); 

(North,  South) 

Barrow  (1983);  Bennett  & 
Tinsley  (1995) 

1957  Daly  City,  California 

5.3 

Marina  District  (2,  3,  4,  5, 

Kayen  et  al.  (1990); 

School) 

Tokimatsu  et  al.  (1991b); 

T.  L.  Youd  (personal 
communication  to  R.  D. 
Andrus,  1999) 

1964  Niigata,  Japan 

7.5 

Niigata  City  (Al,  Cl,  C2, 

Yoshimi  et  al.  (1984; 

Railway  Station) 

1989);  Tokimatsu  et  al. 
(1991a) 

1975  Haicheng,  China 

7.3 

Chemical  Fiber;  Construction 
Building;  Fishery  & 
Shipbuilding;  Glass  Fiber; 
Middle  School;  Paper  Mill 

Arulanandan  et  al.  (1986) 

1979  Imperial  Valley,  Calif. 

6.5 

Heber  Road  (Channel  fill. 

Bennett  et  al.  (1981; 

1981  Westmorland,  Calif. 

5.9 

Point  bar);  Kombloom; 

1984);  Sykora  & Stokoe 

1987  Elmore  Ranch,  Calif 

5.9 

McKim;  Radio  Tower;  Vail 

(1982);  Youd  & Bennett 

1987  Superstition  Hills,  Calif 

6.5 

Canal;  Wildlife 

(1983);  Bierschwale  & 
Stokoe  (1984);  Stokoe  & 
Nazarian  (1984);  Dobry 
et  al.  (1992);  Youd  & 

Holzer  (1994) 

1980  Mid-Chiba,  Japan 

5.9 

Owi  Island  No.  1 

Ishihara  et  al.  (1981;  1987) 

1985  Chiba- Ibaragi-Kenkyo,  Japan 

6.0 

1983  Borah  Peak,  Idaho 

6.9 

Andersen  Bar;  Goddard 

Youd  et  al.  (1985);  Stokoe 

Ranch;  Mackay  Dam 

et  al.  (1988a);  Andrus  et  al. 

Downstream  Toe;  North 

Gravel  Bar;  Pence  Ranch 

(1992);  Andrus  (1994) 

1986  Event  LSST2,  Taiwan 

5.3 

Lotung  LSST  Facility 

Shen  et  al.  (1991); 

Event  LSST3,  Taiwan 

5.5 

EPRI(1992) 

Event  LSST4,  Taiwan 

6.6 

Event  LSST6,  Taiwan 

5.4 

Event  LSST7,  Taiwan 

6.6 

Event  LSST8,  Taiwan 

6.2 

Event  LSST12,  Taiwan 

6.2 

Event  LSST13,  Taiwan 

6.2 

Event  LSST16,  Taiwan 

7.6 

1987  Chiba-Toho-Oki,  Japan 

6.5 

Sunamachi 

Ishihara  et  al.  (1989) 
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Table  3.1  (cont.)  - Earthquakes  and  Sites  Used  to  Establish  Liquefaction  Resistance  Curves. 


Earthquake 

Moment 

Magnitude 

Site 

Reference 

(1) 

(2) 

(3) 

(4) 

1989  Loma  Prieta,  Calif. 

7.0 

Bay  Bridge  Toll  Plaza,  Bay 
Farm  Island  (Dike,  South 

Loop  Road);  Port  of  Oakland; 
Port  of  Richmond 

Stokoe  et  al.  (1992); 

Mitchell  et  al.  (1994) 

Coyote  Creek;  Salinas  River 
(North,  South); 

Barrow  (1983); 

M.  J.  Bennett  (personal 
communication  to  R.  D. 
Andrus,  1995);  Bennett 
and  Tinsley  (1995) 

Marina  District  (2,  3,  4,  5, 
school) 

Kayen  et  al.  (1990); 
Tokimatsu  et  al.  (1991b) 

Moss  Landing  (Harbor 

Office,  Sandholdt  Road, 

State  Beach) 

Boulanger  et  al.  (1995); 
Boulanger  et  al.  (1997) 

Santa  Cruz  (SC02,  SC03, 
SC04,  SC05,  SC13,  SC14) 

Hryciw  (1991); 

Hryciw  et  al.  (1998) 

Treasure  Island  Fire  Station 

Hryciw  et  al.  (1991); 

Redpath  (1991);  Gibbs 
et  al.  (1992);  Furhriman 
(1993);  Andrus  (1994); 
de  Alba  et  al.  (1994) 

Treasure  Island  Perimeter 
(Approach  to  Pier,  UM03, 
UM05,  UM06,  UM09) 

Geomatrix  Consultants 
(1990);  Hryciw  (1991); 

R.  D.  Hryciw  (personal 
communication  to  R.  D. 
Andrus,  1998);  Hryciw 
et  al.  (1998);  Andrus  et  al. 
(1998a,  1998b) 

1993  Kushiro-Oki,  Japan 

8.3 

Kushiro  Port  (2,  D) 

lai  etal.  (1995);  S.  lai 
(personal  communication 
to  R.  D.  Andrus,  1997) 

1993  Hokkaido-Nansei-Oki,  Japan 

8.3 

Pension  House 

Kokusho  et  al.  (1995a, 

1995b, 1995c) 

Hakodate  Port 

S.  lai  (personal 
communication  to 

R.  D.  Andrus,  1997) 

1994  Northridge,  Calif. 

6.7 

Rory  Lane 

Abdel-Haq  & Hryciw 
(1998) 
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Table  3.1  (cont.)  - Earthquakes  and  Sites  Used  to  Establish  Liquefaction  Resistance  Curves. 


Earthquake 

Moment 

Magnitude 

Site 

Reference 

(1) 

(2) 

(3) 

(4) 

1995  Hyogo-Ken  Nanbu,  Japan 

6.9 

Hanshin  Expressway  5 
(3,  10,  14,  25,  29);  Kobe- 
Nishinomiya  Expressway 
(3,  17,  23,  28) 

Hamada  et  al.  (1995); 
Hanshin  Expressway 

Public  Corporation  (1998) 

KNK;  Port  Island  (Downhole 
Array);  SGK 

Sato  et  al.  (1996); 

Shibata  et  al.  (1996) 

Port  Island  (Common 

Factory) 

Ishihara  et  al.  (1997); 
Ishihara  et  al.  (1998) 

Kobe  Port  (7C);  Port  Island 
(1C,  2C) 

Inatomi  et  al.  (1997); 
Hamada  et  al.  (1995) 

Kobe  Port  (LPG  Tank  Yard) 

S.  Yasuda  (personal 
communication  to 

R.  D.  Andrus,  1997) 
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Fig.  3.1  - Relationship  Between  Moment  Magnitude  and  Various  Magnitude  Scales 
(after  Heaton  et  al.,  1982). 
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Values  of  reported  by  the  investigator(s)  are  used  directly.  The  one  exception  is  for 
the  downhole  array  located  at  the  Marina  District  School  site  in  San  Francisco,  California.  A 
reevaluation  of  the  field  data  indicates  that  values  reported  for  the  critical  layer  at  this  site 
are  too  high.  They  are  recalculated  using  the  pseudo-interval  method,  as  discussed  in  Section 
3.2.2. 


Only  the  crosshole  measurements  made  with  shear  waves  having  particle  motion  in  the 
vertical  direction  are  used.  Crosshole  measurements  near  the  critical  layer  boundary  that  seem 
high,  and  could  represent  refracted  waves,  are  not  included  in  the  average. 

Some  V^-values  are  from  measurements  made  before  the  earthquake,  others  following 
the  earthquake.  No  adjustments  are  made  to  compensate  for  changes  in  soil  density  and  due 
to  ground  shaking. 

3.1.3  Measurement  Depth 

In  situ  measurements  may  be  reported  at  discrete  depths  or  for  continuous  intervals, 
depending  on  the  test  method.  When  velocities  are  reported  for  continuous  intervals,  as  is  the 
case  for  downhole  and  SASW  measurements,  the  depth  to  the  center  of  each  interval  is 
assumed.  Thus,  if  the  reported  profile  has  ten  velocity  layers,  it  is  assumed  that  the  profile 
consists  of  ten  “measurements”  with  depths  at  the  center  of  each  layer. 

3.1.4  Case  History 

In  this  report,  a case  history  is  defined  as  a seismic  event  and  a test  array.  For  example, 
at  the  Treasure  Island  Fire  Station  site,  crosshole  measurements  were  made  between  five 
different  pairs  of  boreholes,  downhole  measurements  were  made  by  two  different  investigators, 
seismic  cone  measurements  were  made  at  one  location,  and  SASW  measurements  were  made 
along  one  alignment.  Thus,  a total  of  nine  case  histories  are  identified  for  the  Fire  Station  site 
and  the  1989  Loma  Prieta,  California  earthquake.  At  the  Marina  District  School  site,  downhole 
measurements  were  made  at  one  location.  Estimates  of  ground  surface  acceleration  at  this  site 
are  available  for  the  1957  Daly  City  and  1989  Loma  Prieta  earthquakes.  Thus,  two  case 
histories  are  identified  for  the  Marina  District  School  site.  Combining  the  26  seismic  events 
and  the  139  test  arrays,  a total  of  225  case  histories  are  obtained  with  149  from  the  United 
States,  36  from  Taiwan,  34  from  Japan,  and  6 from  China. 
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The  two  exceptions  to  this  definition  are  the  Owi  Island  No.  1 and  Moss  Landing 
Sandholdt  Road  UC-4  sites  where  additional  subsurface  information  is  available.  At  Owi 
Island,  pore  pressure  transducers  recorded  pore-water  pressure  buildup  for  two  separate  layers. 
At  Moss  Landing,  inclinometer  measurements  indicated  lateral  movement  in  an  upper  loose 
layer  and  no  lateral  movement  in  a lower  dense  layer.  Thus,  two  case  histories  are  identified 
for  each  of  these  two  test  arrays. 

3.1.5  Liquefaction  Occurrence 

It  is  important  to  realize  that  the  occurrence  of  liquefaction,  in  this  evaluation,  is  based 
on  the  appearance  of  surface  evidence,  such  as  sand  boils,  ground  cracks  and  fissures,  and 
ground  settlement.  Case  histories  are  classified  as  non-liquefaction  when  no  liquefaction 
effects  were  observed.  At  the  Owi  Island  No.  1,  Lotung  LSST  Facility,  Sunamachi,  Wildlife 
(1987  earthquakes),  and  Port  Island  sites,  the  assessment  of  liquefaction  or  non-liquefaction 
occurrence  is  supported  by  pore- water  pressure  measurements.  Figure  3.2  shows  the 
distribution  of  case  histories  with  earthquake  magnitude.  Of  the  225  case  histories,  90  are 
liquefaction  case  histories  and  135  are  non-liquefaction  case  histories. 
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Fig.  3.2  - Distribution  of  Liquefaction  and  Non-Liquefaction  Case  Histories  by  Earthquake 
Magnitude. 
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3.1.6  Critical  Layer 


The  critical  layer  is  the  layer  of  non-plastic  soil  below  the  ground  water  table  where 
values  of  as  defined  in  Chapter  2,  and  penetration  resistance  are  generally  the  least,  and 
where  the  cyclic  stress  ratio  relative  to  is  the  greatest.  Figure  3.3  presents  the  cumulative 
relative  frequency  distribution  for  the  case  histories  by  critical  layer  thickness.  Critical  layer 
thicknesses  range  from  1 m to  as  much  as  13  m.  About  50  % of  the  case  histories  have  a 
critical  layer  thickness  less  than  3.5  m;  90  % of  the  case  histories  have  a critical  layer  thickness 
less  than  7 m. 

Figure  3.4  presents  the  cumulative  relative  frequency  distribution  for  the  case  histories 
by  average  measurement  depth  in  the  critical  layer.  The  average  depths  of  the 
measurements  are  between  2 m and  11  m for  nearly  all  case  histories.  Over  50  % of  the  case 
histories  have  average  measurement  depths  less  than  5.5  m.  About  90  % of  the  case  histories 
have  average  measurement  depths  less  than  8 m. 

Materials  comprising  the  critical  layers  range  from  clean  fine  sand  to  sandy  gravel  with 
cobbles  to  profiles  including  silty  clay  layers.  Figure  3.5  summarizes  the  average  fines  content 
(silt  and  clay)  for  the  case  histories  grouped  according  to  earthquake  moment  magnitude.  Of 
the  225  case  histories,  57  are  for  soils  with  5 % or  less  fines,  98  for  soils  with  6 % to  34  % 
fines,  and  70  for  soils  with  35  % or  more  fines.  About  20  % of  the  case  histories  are  for  soils 
containing  more  than  10  % gravel. 

About  70  % of  the  case  histories  are  for  natural  soils  deposits,  with  many  formed  by 
alluvial  processes.  The  other  30  % are  for  hydraulic  or  dumped  fills.  Eight  of  the  fills  have 
been  densified  by  soil  improvement  techniques. 

At  least  85  % of  the  case  histories  are  of  Holocene  age  (<  10  000  years).  While  the  age 
of  the  other  15  % is  unknown,  they  are  believed  to  be  also  of  Holocene  age. 

3.1.7  Ground  Water  Table 

The  ground  water  table  for  nearly  all  case  histories  lies  between  depths  of  0.5  m and 
6 m,  as  shown  in  Fig.  3.6.  Nearly  60  % of  the  case  histories  have  water  table  depths  less  than 
2 m.  About  90  % of  the  case  histories  have  water  table  depths  less  than  4.5  m. 

Artesian  pressures  are  reported  for  the  Lotung  Large-Scale  Seismic  Test  (LSST) 
Facility  site  in  Taiwan.  At  this  site,  the  pore-water  pressure  distribution  is  assumed  to  vary 
linearly  from  a pressure  head  of  8.1  m at  a depth  of  7 m to  a pressure  head  of  1.9  m at  a depth 
of  2 m. 
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Fig.  3.3  - Cumulative  Relative  Frequency  of  Case  History  Data  by  Critical  Layer  Thickness. 
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Fig.  3.4  - Cumulative  Relative  Frequency  of  Case  History  Data  by  Average  Depth  of 
Measurements  in  Critical  Layer. 
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3.1.8  Total  and  Effective  Overburden  Stresses 


Values  of  total  and  effective  overburden  stresses  are  estimated  using  densities  reported 
by  the  investigator(s).  When  no  densities  are  reported,  typical  values  for  soils  with  similar 
grain  size,  penetration  and  velocity  characteristics  are  assumed.  In  most  instances,  the  assumed 
densities  are  1.76  Mg/m^  for  soils  above  the  water  table  and  1.92  Mg/m^  for  soils  below  the 
water  table. 

3.1.9  Average  Peak  Ground  Acceleration 

Average  values  of  peak  horizontal  ground  surface  acceleration,  a^,  are  determined  by 
averaging  estimates  reported  by  the  investigator(s)  and  estimates  made  as  part  of  this  study 
using  attenuation  relationships  developed  from  published  ground  surface  acceleration  data. 
Since  many  published  attenuation  relationships  are  based  on  both  peak  values  obtained  from 
ground  motion  records  for  the  two  horizontal  directions  (sometimes  referred  to  as  the  randomly 
oriented  horizontal  component),  the  geometric  mean  (square  root  of  the  product)  of  the  two 
peak  values  is  used.  Use  of  the  geometric  mean  is  consistent  with  the  development  of  the  SPT- 
based  procedure  (Youd  et  al.,  1997).  For  the  cases  in  this  study,  the  difference  between  the 
geometric  mean  and  arithmetic  mean  values  is  generally  small,  within  about  5 %. 

3.1.10  Average  Cyclic  Stress  Ratio 

Cyclic  stress  ratios,  CSR,  are  first  calculated  for  each  “measurement”  depth  within  the 
critical  layer  using  Eq.  2.1,  and  then  averaged.  Values  of  are  estimated  using  the  average  by 
Seed  and  Idriss  (1971)  shown  in  Fig.  2.1.  These  values  are  used  to  follow  the  traditional 
format  of  the  SPT-  and  CPT-based  procedures  where  the  magnitude  scaling  factor  is  used  to 
account  for  all  effects  of  earthquake  magnitude. 

3.1.11  Average  Overburden  Stress-Corrected  Shear  Wave  Velocity 

Values  of  Vg  within  the  critical  layer  are  first  corrected  for  overburden  stress  using  Eq. 
1 .4,  and  then  averaged.  The  number  of  values  included  in  the  average  range  from  1 to  22  (see 
Appendix  C).  About  80  % of  the  case  histories  have  2 to  7 values  in  the  average.  No 
adjustments  are  made  for  possible  variations  between  seismic  test  methods  due  to  different 
source-receiver  orientations  with  respect  to  the  stress  state  in  the  soil.  In  the  calculations,  each 
site  is  assumed  to  be  level  ground. 


41 


3.2  SAMPLE  CALCULATIONS 


Calculations  for  two  sites  shaken  by  the  1989  Loma  Prieta,  California  earthquake  (M^  = 
7.0)  are  presented  below  to  illustrate  how  values  of  CSR  and  overburden  stress-corrected  shear 
wave  velocity,  are  determined.  The  two  sites  are  Treasure  Island  Fire  Station  and  Marina 
District  School. 

3.2.1  Treasure  Island  Fire  Station 

Treasure  Island  is  a man-made  island  located  in  the  San  Francisco  Bay  along  the  Bay 
Bridge  between  the  cities  of  San  Francisco  and  Oakland.  It  was  constructed  in  1936-37  by 
hydraulic  filling  behind  a perimeter  rock  dike. 

Extensive  field  tests  have  been  conducted  at  the  fire  station  on  Treasure  Island.  Figure 
3.7  presents  two  profiles  for  the  site.  The  profile  determined  by  crosshole  testing  is  from 
Fuhriman  (1993).  The  other  profile  is  based  on  unpublished  SASW  test  results  by  The 
University  of  Texas  at  Austin  in  1992.  Also  presented  in  Fig.  3.7  is  the  soil  profile  for  the  site. 
From  the  description  by  de  Alba  et  al.  (1994),  the  upper  4.5  m of  soil  consists  of  silty  sand  fill, 
possibly  formed  by  dumping.  Between  depths  of  4.5  m and  12.2  m,  the  soil  consists  of  silty 
sand  to  clayey  sand,  formed  by  hydraulic  filling.  Beneath  the  hydraulic  fill  are  natural  clayey 
soils.  The  ground  water  table  lies  near  the  ground  surface  at  a depth  of  1.4  m.  The  critical 
layer  is  determined  to  be  between  depths  of  4.5  m and  7 m,  where  the  soil  is  non-plastic,  lies 
below  the  water  table,  and  exhibits  the  lowest  values  of  relative  to  the  highest  values  of 
CSR  in  the  layer  (see  Fig.  5.1). 

During  the  1989  Loma  Prieta  earthquake,  a seismograph  station  at  the  fire  station 
recorded  ground  surface  accelerations.  The  peak  values  in  the  two  horizontal  accelerometer 
records  are  0.16  g and  0.11  g (Brady  and  Shakal,  1994).  Unlike  recordings  at  other 
seismograph  stations  located  on  soft-soils  in  the  Bay  area,  there  is  a sudden  drop  in  the 
recorded  acceleration  at  about  15  seconds  and  small  motion  afterwards  (Idriss,  1990).  De  Alba 
et  al.  (1994)  attribute  this  behavior  to  liquefaction  of  an  underlying  sand.  However,  no  sand 
boils  or  ground  cracks  occurred  at  the  site.  The  nearest  liquefaction  effect  observed  is  a sand 
boil  located  100  m from  the  site  (Geometric  Consultants,  1990;  Bennett,  1994;  Power  et  al., 
1998).  Thus,  this  site  is  classified  as  a non-liquefaction  site  during  this  event  by  the  definition 
given  in  Section  3.1.4. 
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Fig.  3.7  - Shear  Wave  Velocity  and  Soil  Profiles  for  the  Treasure  Island  Fire  Station  Site. 
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Sample  calculations  for  the  crosshole  and  SASW  test  arrays  are  summarized  in  Tables 
3.2  and  3.3,  respectively.  The  data  points  used  in  the  calculations  are  shown  by  the  open 
symbols  in  Fig.  3.7.  Total  and  effective  overburden  stresses  are  calculated  assuming  densities 
of  1.76  Mg/m^  above  the  water  table  and  1.92  Mg/m^  below  the  water  table.  Stress  reduction 
coefficients  are  estimated  using  the  average  curve  by  Seed  and  Idriss  (1971)  shown  in  Fig.  2.1. 
The  geometric  mean  of  the  two  peak  values  observed  in  the  horizontal  ground  surface 
acceleration  records  is  0.13  g.  Using  these  parameters,  values  of  CSR  and  are  calculated  for 
the  crosshole  measurement  at  depth  of  4.6  m as  follows: 


CSR 


= 0.65 
= 0.65 


f \ 

Clmax 

^(7v^ 

l ^ J 

Vcr'vj 

"0.13gY84.0^ 

V 5 JI52.7J 


0.97  =0.131 


(3,1) 


and 


SI 


K(y'vJ 


= 134 


.52.7J 


0.25 


= 158  m/s 


(3.2) 


Representative  values  of  CSR  and  used  to  defined  the  two  case  histories  are  determined  by 
averaging  values  for  each  “measurement”  depth  within  the  critical  layer,  as  shown  in  Tables  3.2 
and  3.3. 


3.2.2  Marina  District  School 


Downhole  seismic  tests  were  conducted  at  the  Winfield  Scott  School  in  the  Marina 
District  of  San  Francisco  by  Kayen  et  al.  (1990).  Figure  3.8  presents  soil  and  velocity  profiles 
for  the  site.  The  critical  layer  lies  between  depths  of  2.7  m,  the  ground  water  table  depth,  and 
4.3  m,  the  base  of  sand  fill.  The  average  profile  shown  in  Fig.  3.8  was  determined  by  Kayen 
et  al.,  and  was  based  on  best-fit  line  segments  through  travel  time  measurements  plotted  versus 
depth.  The  second  profile  is  determined  using  the  pseudo-interval  method  (This  Report),  as 
illustrated  in  Fig.  3.9.  Both  methods  should  provide  similar  average  values  over  the  same  depth 
interval.  However,  the  layering  assumed  for  the  best-fit  line  segment  method  does  not  seem 
appropriate  for  the  fill.  For  this  reason,  values  of  based  on  the  pseudo-interval  method  are 
used  in  this  analysis. 
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Table  3.2  - Sample  Calculations  for  the  Treasure  Island  Fire  Station  Site,  Crosshole  Test  Array 
Bl-  B4,  and  the  1989  Loma  Prieta  Earthquake. 


Measurement 

Number 

(1) 

Average 

Depth, 

m 

(2) 

Measured 

Shear 

Wave 

Velocity, 

m/s 

(3) 

Total 

Overburden 

Stress', 

kPa 

(4) 

Effective 

Overburden 

Stress', 

kPa 

(5) 

Stress 

Reduction 

Coefficient, 

rd 

(6) 

Cyclic 

Stress 

Ratio^ 

CSR 

(7) 

Overburden 

Stress- 

Corrected 

Shear  Wave 
Velocity, 
m/s 

(8) 

1 

4.57 

134 

84.0 

52.7 

0.97 

0.13 

158 

2 

5.49 

133 

111.3 

60.9 

0.96 

0.14 

150 

3 

6.40 

144 

118.5 

69.2 

0.95 

0.14 

158 

Average 

5.5 

137 

101.3 

60.9 

0.96 

0.14 

155 

'Assuming  water  table  at  1.4  m;  and  material  densities  are  1.76  Mg/m^  above  the  water  table 


and  1 .92  Mg/m^  below  the  water  table. 

^Based  on  average  values  determined  by  Seed  and  Idriss  (1971). 
Assuming  peak  horizontal  ground  surface  acceleration  is  0.13  g. 


Table  3.3  - Sample  Calculations  for  the  Treasure  Island  Fire  Station  Site,  SASW  Test  Array, 
and  the  1989  Loma  Prieta  Earthquake. 


Measured 

Overburden 

Shear 

Stress 

Cyclic 

Stress- 

Wave 

Total 

Effective 

Reduction 

Stress 

Corrected 

Measurement 

Average 

Velocity, 

Overburden 

Overburden 

Coefficient, 

Ratio^ 

Shear  Wave 

Number 

Depth, 

F„ 

Stress', 

Stress', 

O 

CSR 

Velocity,  V^j, 

m 

m/s 

kPa 

kPa 

m/s 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

1 

5.34 

131 

98.4 

59.6 

0.96 

0.14 

149 

Average 

5.3 

131 

98.4 

59.6 

0.96 

0.14 

149 

'Assuming  water  table  at  1.4  m;  and  material  densities  are  1.76  Mg/m^  above  the  water  table 


and  1.92  Mg/m^  below  the  water  table. 

^Based  on  average  values  determined  by  Seed  and  Idriss  (1971). 
Assuming  peak  horizontal  ground  surface  acceleration  is  0.13  g. 
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Fig.  3.8  - Shear  Wave  Velocity  and  Soil  Profiles  for  the  Marina  District  School  Site 
(Kayen  et  al.,  1990). 


46 


Fig.  3.9  - General  Configuration  of  the  Downhole  Seismic  Test  Using  the  Pseudo-Interval 
Method  to  Calculate  Shear  Wave  Velocity. 


Many  structures,  pavements,  and  public  works  near  the  school  sustained  heavy  damage 
during  the  1989  earthquake  (Kayen  et  al.,  1990).  This  damage  was  due  to  liquefaction  of  the 
sand  fill.  From  maps  prepared  by  Pease  and  O’Rourke  (1995),  the  site  lies  on  the  margin  of  the 
1906  water  front  and  artificial  fill  where  about  40  mm  of  settlement  occurred.  Mapped  sand 
boils  and  ground  cracks  lie  just  east  of  the  site.  Based  on  these  observations,  this  site  is 
classified  as  a liquefaction  site  during  this  earthquake. 

The  Marina  District  and  Treasure  Island  are  located  about  82  km  from  the  1989  surface 
rupture.  The  geometric  mean  value  of  recorded  at  Treasure  Island  during  the  earthquake  is 
0.13  g.  The  attenuation  relationship  by  Idriss  (1991)  for  1989  strong  ground  motion  records 
from  soft-soil  sites  provides  a median  value  of  0.16  g,  assuming  a distance  of  82  km  from 
surface  rupture.  Thus,  a peak  horizontal  ground  surface  acceleration  of  0.15  g,  the  average  of 
these  two  estimates,  is  assumed  in  the  analysis. 

Sample  calculations  for  the  Marina  District  School  site  are  summarized  in  Table  3.4. 
The  locations  of  measurements  are  assumed  midway  between  receiver  positions,  as  shown  in 
Fig.  3.8.  Total  and  effective  overburden  stresses  are  estimated  assuming  densities  of 
1.76  Mg/m^  above  the  water  table  and  1.92  Mg/m^  below  the  water  table.  The  ground  water 
table  is  at  a depth  of  about  2.7  m.  Average  values  of  CSR  and  defining  the  case  history  are 
determined  by  averaging  values  for  the  two  “measurement”  depths,  as  shown  in  Table  3.4 
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3.3  SUMMARY 


The  case  history  data  described  in  this  chapter  are  limited  to  level  and  gently  sloping 
sites  with  the  following  characteristics: 

(1)  average  critical  layer  depths  less  than  10  m; 

(2)  uncemented  soils  of  Holocene  age; 

(3)  ground  water  table  depths  between  0.5  m and  6 m;  and 

(4)  all  V5  measurements  from  below  the  water  table. 

Of  the  225  case  histories,  57  are  for  soils  with  5 % or  less  fines  (silt  and  clay),  98  for  soils  with 
6 % to  34  % fines,  and  70  for  soils  with  35  % or  more  fines.  About  20  % of  the  case  histories 
are  for  soils  containing  more  than  10  % gravel.  Nearly  50  % of  the  case  histories  are  for 
earthquake  magnitudes  near  7. 


Table  3.4  - Sample  Calculations  for  the  Marina  District  School  Site  and  the  1989  Loma  Prieta 
Earthquake. 


Measurement 

Number 

(1) 

Average 

Depth, 

m 

(2) 

Measured 

Shear 

Wave 

Velocity', 

Vs, 

m/s 

(3) 

Total 

Overburden 

Stress^, 

kPa 

(4) 

Effective 

Overburden 

Stress^ 

kPa 

(5) 

Stress 

Reduction 

Coefficient, 

0 

(6) 

Cyclic 

Stress 

Ratio", 

CSR 

(7) 

Overburden 

Stress- 

Corrected 

Shear  Wave 

Velocity,  V^y, 

m/s 

(8) 

1 

3.02 

87 

52.6 

49.9 

0.98 

0.10 

104 

2 

3.94 

136 

70.0 

58.2 

0.97 

0.11 

156 

Average 

3.5 

112 

61.3 

54.1 

0.98 

0.11 

130 

'Based  on  pseudo-interval  method. 

^Assuming  water  table  at  2.7  m;  and  material  densities  are  1.76  Mg/m^  above  the  water  table 
and  1.92  Mg/m^  below  the  water  table. 

^Based  on  average  values  determined  by  Seed  and  Idriss  (1971). 

'‘Assuming  peak  horizontal  ground  surface  acceleration  is  0.15  g. 
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CHAPTER  4 


LIQUEFACTION  RESISTANCE  FROM  CASE  HISTORY  DATA 


To  establish  the  recommended  curves  for  liquefaction  resistance  evaluation,  cyclic 
stress  ratios  and  shear  wave  velocities  for  the  225  case  histories  described  in  Chapter  3 are 
plotted  in  Figs.  4.1  through  4.6  for  earthquakes  with  moment  magnitudes  ranging  from  5.5 
through  8,  respectively.  The  plotted  data  have  been  separated  into  three  categories:  (1)  sands 
and  gravels  with  average  fines  (silt  and  clay)  content  less  than  or  equal  to  5 %;  (2)  sand  and 
gravels  with  average  fines  content  of  6 % to  34  %;  and  (3)  sands  and  silts  with  average  fines 
content  greater  than  or  equal  to  35  %.  Also  shown  in  Figs.  4.1  through  4.6  are  the 
recommended  liquefaction  resistance  curves  based  on  Eq.  2.21.  This  chapter  presents  the 
development  of  these  curves. 


4.1  LIMITING  UPPER  V^j  VALUE  FOR  LIQUEFACTION  OCCURRENCE 

Figures  4.1  through  4.6  show  that  the  available  case  history  data  above  a cyclic  stress 
ratio  of  about  0.35  are  limited.  Thus,  current  estimates  of  Vs\  rely  in  part  on  penetration-shear 
wave  velocity  correlations. 

In  the  SPT-based  procedure,  a corrected  blow  count  of  30  is  assumed  as  the  limiting 
upper  value  for  liquefaction  occurrence  in  clean  sands  above  CRR  of  0.6  and  magnitude  7.5 
earthquakes  (Seed  et  al.,  1985;  Youd  et  al.,  1997).  Table  4.1  presents  estimates  of  equivalent 
V^j  for  corrected  blow  count  of  30.  The  relationship  by  Ohta  and  Goto  (1978)  modified  to  blow 
count  with  theoretical  free-fall  energy  of  60  % (Seed  et  al.,  1985)  suggests  equivalent  V^j 
values  of  207  m/s  for  Holocene  sands  and  227  m/s  for  Holocene  gravels,  assuming  that  a depth 
of  10  m is  equivalent  to  an  overburden  stress  of  100  kPa.  The  stress-corrected  crosshole 
measurements  collected  by  Sykora  (1987b)  for  Holocene  sands  and  non-plastic  silty  sands 
below  the  ground  water  table  with  corrected  blow  count  between  25  and  35  exhibit  an  average 
value  of  206  m/s  and  a standard  deviation  of  41  m/s.  The  correlation  by  Rollins  et  al.  (1998a) 
provides  a best-fit  value  of  232  m/s  for  Holocene  gravels  with  corrected  blow  count  of  30. 
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Fig.  4.1  - Case  History  Data  for  Earthquakes  with  Magnitude  Near  5.5  Based  on 

Overburden  Stress-Corrected  Shear  Wave  Velocity  and  Cyclic  Stress  Ratio 
with  Recommended  Liquefaction  Resistance  Curves. 
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Fig.  4.2  - Case  History  Data  for  Earthquakes  with  Magnitude  Near  6 Based  on 

Overburden  Stress-Corrected  Shear  Wave  Velocity  and  Cyclic  Stress  Ratio 
with  Recommended  Liquefaction  Resistance  Curves. 
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Fig.  4.3  - Case  History  Data  for  Earthquakes  with  Magnitude  Near  6.5  Based  on 

Overburden  Stress-Corrected  Shear  Wave  Velocity  and  Cyclic  Stress  Ratio 
with  Recommended  Liquefaction  Resistance  Curves. 
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Fig.  4.4  - Case  History  Data  For  Earthquakes  with  Magnitude  Near  7 Based  on 

Overburden  Stress-Corrected  Shear  Wave  Velocity  and  Cyclic  Stress  Ratio 
with  Recommended  Liquefaction  Resistance  Curves. 
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Fig.  4.5  - Case  History  Data  for  Earthquakes  with  Magnitude  Near  7.5  Based  on 

Overburden  Stress-Corrected  Shear  Wave  Velocity  and  Cyclic  Stress  Ratio 
with  Recommended  Liquefaction  Resistance  Curves. 
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Fig.  4.6  - Case  History  Data  for  Earthquakes  with  Magnitude  Near  8 Based  on 

Overburden  Stress-Corrected  Shear  Wave  Velocity  and  Cyclic  Stress  Ratio 
with  Recommended  Liquefaction  Resistance  Curves. 
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Table  4.1  - Estimates  of  Equivalent  for  Holocene  Sands  and  Gravels  Below  the  Ground 
Water  Table  with  Corrected  SPT  Blow  Count  of  30. 


Reference 

(1) 

Relationship 

(2) 

Equivalent 

Vo;  Estimate 
(m/s) 

(3) 

Assumptions 

(4) 

Ohta  & 

Goto 

(1978);  also 
given  in 
report  by 
Sykora 
(1987a, 
page  29) 

Vj  = 69  F,F2 

Nj  = SPT  blow  count 

measured  in  Japanese 
practice 
z = depth,  m 

Fj  = 1 .00  for  Holocene- 
age  soils 

F2  = 1.085  for  sands; 

1.189  for  gravel 

...best-fit  relationship  for  289 
sets  of  SPT  and 
measurements  from  Japan 

207 

...for  Holocene 
sands 

227 

...for  Holocene 
gravels 

1.  Nj  = 610161 

2.  a4,=  30 

3.  z = 10  m is  equivalent  to 
an  overburden  stress  of 

100  kPa 

4.  All  measurements  are 
from  below  the  ground 
water  table 

Sykora 
(1987b, 
page  90); 
This  Report 

Correlation  between  {Nj)^q  and 
crosshole  V^,  normalized  to 
effective  overburden  stress, 
measurements  for  Holocene 
sands  and  non-plastic  silty 
sands  below  the  ground  water 
table  at  sites  in  U.S.A.;  16  sets 
of  measurements  (with  known 
SPT  equipment) 

206 

...for  Holocene 
sands  and  non- 
plastic silty 
sands  below  the 
water  table 
...standard 
deviation  is 

41  m/s 

1 . Average  for  V^j  values 
with  {Nj)^q  between 

25  and  35 

2.  o'v  = 100  kPa 

Rollins  et 
al.  (1998a) 

V,  = 53  (AT,/  {o'T" 

...best-fit  relationship  using 
equivalent  values  from 

Becker  Penetration  Tests  and 
measurements;  186  points 
from  7 Holocene  gravel  sites 

232 

...for  Holocene 
gravels 
...most  of  data 
lie  within 
±25  % of 
relationship 

1.  N,,  = 00 

2.  cr'v  = 100  kPa 

3.  All  measurements  are 
from  below  the  ground 
water  table 

This  Report 
(see  Fig. 

4.7) 

ysi-B,  [{N,)X 

Bj  =93.2  ±6.5 

B2  =0.231  ±0.022 

...best-fit  relationship  for 
uncemented,  Holocene-age 
sands  with  less  than  10  % non- 
plastic fines;  25  sets  of  average 
SPT  and  measurements  all 
from  below  the  water  table 

204 

...for  Holocene 
clean  sands 
below  the  water 
table 
...residual 
standard 
deviation  is 

12  m/s 

1 . Average  for 

with  = 30 

2.  a'v  = 100  kPa 

3.  Corrected  blow  count 
based  on  procedures 
given  in  Seed  et  al. 

(1985)  and  Robertson 
and  Wride  (1997;  1998) 
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Rollins  et  al.  noted  that  the  majority  of  the  data  analyzed  fall  within  ±25  % of  their  correlation. 
Figure  4.7  presents  average  values  of  and  corrected  blow  count  for  soil  layers  with  less  than 
10  % fines  at  several  sites  listed  in  Table  3.1.  Also  shown  in  Fig.  4.7  is  the  best-fit  relationship 
for  the  plotted  data.  The  plotted  data  exhibit  a mean  value  of  204  m/s  at  a corrected  blow 
count  of  30  and  a residual  standard  deviation,  of  12  m/s. 

In  the  CPT-based  procedure,  a normalized  cone  tip  resistance  of  160  is  assumed  as  the 
limiting  upper  value  for  liquefaction  occurrence  in  clean  sand  above  CRR  of  0.6  and  magnitude 
7.5  earthquakes  (Youd  et  al.,  1997;  Robertson  and  Wride,  1998).  Figure  4.8  presents  average 
values  of  and  normalized  tip  resistance  for  soil  layers  with  less  than  10  % fines  at  several 
sites  listed  in  Table  3.1.  Also  shown  in  Fig.  4.8  is  the  best-fit  relationship  for  the  plotted  data. 
Table  4.2  summarizes  the  general  characteristics  of  the  plotted  data  and  best-fit  relationship. 
The  plotted  data  exhibit  a mean  value  of  193  m/s  at  a normalized  tip  resistance  of  160  and  a 
residual  standard  deviation  of  19  m/s. 

From  these  estimates,  a value  of  210  m/s  is  assumed  equivalent  to  a corrected  blow 
count  of  30  and  normalized  tip  resistance  of  160  in  clean  sands.  A limiting  upper  value  of 
210  m/s  for  liquefaction  occurrence  is  less  than  the  general  consensus  value  of  230  m/s  agreed 
upon  at  the  1998  MCEER  workshop.  Figures  4.7  and  4.8  have  been  added  to  this  report  to 
provide  additional  evidence  to  support  the  use  of  210  m/s  as  the  limiting  upper  value  of  for 
clean  sands  and  gravels. 


Table  4.2  - Estimates  of  Equivalent  for  Holocene  Sands  and  Gravels  Below  the  Ground 
Water  Table  with  Normalized  Cone  Tip  Resistance  of  160. 


Reference 

(1) 

Relationship 

(2) 

Equivalent 

Vo,  Estimate 
(m/s) 

(3) 

Assumptions 

(4) 

This  Report 

193 

1 . Average  for  V^j 

(see  Fig. 

...for  Holocene 

with  160 

4.8) 

Bj  = 88.2  ± 15.5 
=0.154  ±0.037 

...best-fit  relationship  for 
uncemented,  Holocene-age 
sands  with  less  than  10  % non- 
plastic fines;  23  sets  of  average 
SPT  and  measurements  all 
from  below  the  water  table 

clean  sands 
below  the  water 
table 
...residual 
standard 
deviation  is 

19  m/s 

2.  cr'v  = 100  kPa 

3.  Normalized  tip 
resistance  based  on 
procedures  given  in 
Robertson  and  Wride 
(1997; 1998) 
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Fig.  4.7  - Variations  in  with  (Nj)^q  for  Uncemented,  Holocene-age  Sands  with  Less 
than  10  % Non-Plastic  Fines. 
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Fig.  4.8  - Variations  in  V^j  with  q^jj^  for  Uncemented,  Holocene-age  Sands  with  Less  than 
10  % Non-Plastic  Fines. 
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For  soils  with  fines  content  greater  than  35  %,  the  SPT-based  chart  by  Seed  et  al.  (1985) 
indicates  a limiting  upper  corrected  blow  count  of  about  2 1 for  liquefaction  occurrence  above 
CRR  of  0.6.  Table  4.3  presents  estimates  of  equivalent  for  blow  count  of  21.  The 
relationship  by  Ohta  and  Goto  (1978)  suggests  equivalent  V^j  values  of  195  m/s  for  Holocene 
sands  and  214  m/s  for  Holocene  gravels.  The  stress-corrected  crosshole  collected  by  Sykora 
(1987b)  for  Holocene  sands  and  non-plastic  silty  sands  below  the  ground  water  table  with 
corrected  blow  count  between  16  and  26  exhibit  an  average  value  of  199  m/s  and  a standard 
deviation  of  36  m/s.  The  correlation  by  Rollins  et  al.  (1998a)  provides  a best-fit  value  of  217 
m/s  for  Holocene  gravels.  From  these  estimates,  a value  of  195  m/s  is  assumed  equivalent 
to  a corrected  blow  count  of  21  in  non-plastic  soils  with  fines  content  > 35  %. 

To  permit  the  liquefaction  resistance  curves  for  magnitude  7.5  earthquakes  shown  in 
Fig.  4.5  to  have  values  between  195  m/s  and  210  m/s  at  CRR  near  0.6,  values  of  Vsi  are 
assumed  to  range  from  200  m/s  to  215  m/s,  respectively.  The  relationship  between  V*si  and 
fines  content  can  be  expressed  by: 

yJi  =215  m/s  for  sands  and  gravels  with  FC  <5  % (4.1a) 

Vsi  = 215  - 0.5(FC-5)  m/s  for  sands  and  gravels  with  5 % < FC  <35  % (4.1b) 

Vsi  = 200  m/s  for  sands  and  silts  with  FC  > 35  % (4.1c) 

where 

FC  = the  average  fines  content  in  percent  by  mass. 

The  CRR-V^j  curves  shown  in  Figs.  4.1  through  4.6  provide  reasonable  bounds  for  the 
case  history  data  shown  above  a cyclic  stress  ratio  of  0.35,  further  supporting  the  use  of  Vsi 
values  given  in  Eq.  4.1. 
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Table  4.3  - Estimates  of  Equivalent  V^j  for  Holocene  Sands  and  Gravels  Below  the  Ground 
Water  Table  with  Corrected  SPT  Blow  Count  of  21. 


Reference 

(1) 

Relationship 

(2) 

Equivalent 
y..  Estimate 
(m/s) 

(3) 

Assumptions 

(4) 

Ohta  & 

Goto 

(1978);  also 
given  in 
report  by 
Sykora 
(1987a, 
page  29) 

= 69  FjF^ 

Nj  = SPT  blow  count 

measured  in  Japanese 
practice 
z = depth,  m 

F;  = 1 .00  for  Holocene- 
age  soils 

F2  = 1.085  for  sands; 

1.189  for  gravel 

195 

...for  Holocene 
sands 

214 

...for  Holocene 
gravels 

1.  A.  = 60/67  A, ^ 

2.  a;,  = 21 

3.  z = 10  m is  equivalent  to 
an  overburden  stress  of 

100  kPa 

4.  All  measurements  are 
from  below  the  ground 
water  table 

...best-fit  relationship  for  289 
sets  of  SPT  and 
measurements  from  Japan 

Sykora 
(1987b, 
page  90); 
This  Report 

Correlation  between  and 

crosshole  V^,  normalized  to 
effective  overburden  stress, 
measurements  for  Holocene 
sands  and  non-plastic  silty 
sands  below  the  water  table  at 
sites  in  U.S.A.;  31  sets  of 
measurements  (with  known 

SPT  equipment) 

199 

...for  Holocene 
sands  and  non- 
plastic silty 
sands  below  the 
water  table 
...standard 
deviation  is 

36  m/s 

1 . Average  for  y^^  values 
with  (A;)gQ  between 

16  and  26 

2.  (j'v  = 100  kPa 

Rollins  et 
al.  (1998a) 

V,  = 53  {a',r 

...best-fit  relationship  using 
equivalent  7/^^-values  from 
Becker  Penetration  Tests  and 
y^  measurements;  186  points 
from  7 Holocene  gravel  sites 

217 

...for  Holocene 
gravels 
...most  of  data 
lie  within 
±25  % of 
relationship 

1.  = 

2.  cj'v  = 100  kPa 

3.  All  measurements  are 
from  below  the  ground 
water  table 
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4.2  CURVE  FITTING  PARAMETERS,  a AND  b 


The  curve  fitting  parameters  a andZ?  in  Eq.  2.21  can  be  approximated  from  the  case 
history  data  assuming  the  values  of  Vsi  given  in  Eq.  4. 1 and  a MSF  relationship.  By  assuming 
a MSF  relationship,  one  can  adjust  the  cyclic  stress  ratios  for  the  case  history  data  to  magnitude 
7.5  as  follows: 


CSR 


1.5  ~ 


CSR 

MSF 


(4.2) 


where 


CSR  = the  cyclic  stress  ratio,  and 

CSR^  5 = the  CSR  corrected  to  magnitude  7.5. 

Three  MSF  relationships  representing  the  range  of  proposed  magnitude  scaling  factors  are 
considered  below  to  establish  the  values  of  a and  b. 

4.2.1  Magnitude  Scaling  Factors  Recommended  by  1996  NCEER  Workshop 

4.2.1.1  Lower  Bound  of  Recommended  Range — Figure  4.9  presents  the  case  history 
data  for  magnitude  5.9  to  8.3  earthquakes  adjusted  using  the  lower  bound  for  the  range  of 
magnitude  scaling  factors  recommended  by  the  1996  NCEER  workshop  (Youd  et  al.,  1997). 
The  lower  bound  is  defined  by  Eq.  2.13  with  n = -2.56,  as  discussed  in  Section  2.4.1.  Also 
shown  in  Fig.  4.9  are  three  liquefaction  resistance  curves  for  earthquakes  with  magnitude  near 
7.5  and  various  fines  content.  The  three  curves  were  determined  through  an  iterative  process  of 
varying  the  values  of  a and  b until  nearly  all  the  liquefaction  case  histories  were  bound  by  the 
curves  with  the  least  amount  of  non-liquefaction  case  histories  in  the  liquefaction  region.  Of 
the  90  liquefaction  case  histories,  only  two  case  histories  incorrectly  lie  in  the  no  liquefaction 
region.  The  final  values  of  a and  b used  to  draw  the  curves  in  Fig.  4.9  are  0.022  and  2.8, 
respectively. 
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Cyclic  Stress  or  Resistance  Ratio,  CSR  or  CRR 
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Fig.  4.9  - Curves  Recommended  for  Calculation  of  CRR  from  Shear  Wave  Velocity 
Measurements  Along  with  Case  History  Data  Based  on  Lower  Bound 
Values  of  MSF  for  the  Range  Recommended  by  the  1996  NCEER 
Workshop  (Youd  et  al.,  1997)  and  Developed  by  Seed  and  Idriss  (1971). 
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The  two  liquefaction  case  histories  in  Fig.  4.9  that  lie  in  the  no  liquefaction  region  are 
for  the  Treasure  Island  UM05  and  UM09  sites.  The  UM05  and  UM09  sites  are  located  along 
the  perimeter  of  Treasure  Island.  Mapped  liquefaction  effects  generated  by  the  1989  Loma 
Prieta  earthquake  near  the  UM05  site  are  ground  cracks  with  50  to  90  mm  of  horizontal 
displacement  (R.  D.  Hryciw,  personal  communication  to  R.  D.  Andrus,  1998;  Power  et  al., 
1998).  The  nearest  mapped  sand  boil  is  located  60  m away  from  the  site.  At  the  UM09  site,  as 
much  as  90  mm  of  vertical  displacement  was  observed  adjacent  to  a building  located  60  m 
inland  from  the  site.  These  displacements  are  small  compared  to  the  meters  of  displacement 
that  are  expected  to  occur  during  larger  ground  shaking.  Thus,  liquefaction  was  marginal  at  the 
UM05  and  UM09  sites,  and  sloping  ground  may  have  been  a factor.  It  should  be  noted  that 
similar  evaluations  are  obtained  using  the  SPT  and  CPT  data  for  these  sites.  The  SPT-  and 
CPT-based  evaluations  for  the  UM05  site  are  discussed  in  Section  4.3.1. 

4.2.1.2  Upper  Bound  of  Recommended  Range — Figure  4.10  presents  the  case  history 
data  for  magnitude  5.9  to  8.3  earthquakes  adjusted  using  the  upper  bound  for  the  range  of 
magnitude  scaling  factors  recommended  by  the  1996  NCEER  workshop  (Youd  et  al.,  1997). 
The  upper  bound  is  defined  by  Eq.  2.13  with  n = -3.3,  as  discussed  in  Section  2.4.1.  Also 
shown  in  Fig.  4.10  are  the  same  three  curves  from  Fig.  4.9.  Many  case  history  data  plot  lower 
in  Fig.  4.10  than  in  Fig.  4.9,  since  the  earthquake  magnitude  is  less  than  7.5  for  most  of  the 
data.  The  downward  shift  in  the  liquefaction  data  points  near  the  curves  at  CRR  of  about  0.08 
is  less  than  0.01.  This  difference  is  not  significant,  and  is  within  the  accuracy  of  the  plotted 
case  history  data. 

4.2.2  Revised  Mlagnitude  Scaling  Factors  Proposed  by  Idriss  (1999) 

Figure  4.1 1 presents  the  case  history  data  for  magnitude  5.9  to  8.3  earthquakes  adjusted 
using  the  revised  magnitude  scaling  factors  and  stress  reduction  coefficients  proposed  by  Idriss 
(1999  as  discussed  in  Section  2.4.2.  Also  shown  in  Fig.  4.11  are  the  same  three  liquefaction 
resistance  curves  from  Fig.  4.9.  Many  of  the  case  history  data  shown  in  Fig.  4.1 1 plot  higher 
than  case  history  data  in  Fig.  4.9,  since  the  earthquake  magnitude  is  less  than  7.5  for  most  of 
the  data.  The  upward  shift  in  the  liquefaction  data  points  near  the  curves  at  CRR  of  about  0.08 
is  less  than  0.01.  This  difference  is  also  not  significant,  and  is  within  the  accuracy  of  the 
plotted  case  history  data. 
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Cyclic  Stress  or  Resistance  Ratio,  CSR  or  CRR 
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Fig.  4.10  - Curves  Recommended  for  Calculation  of  CRR  from  Shear  Wave  Velocity 

Measurements  Along  with  Case  History  Data  Based  on  Upper  Bound  Values 
of  MSF  for  the  Range  Recommended  by  the  1996  NCEER  Workshop  (Youd 
et  al.,  1997)  and  rj  Developed  by  Seed  and  Idriss  (1971). 
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Cyclic  Stress  or  Resistance  Ratio,  CSR  or  CRR 
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Fig.  4.1 1 - Curves  Recommended  for  Calculation  of  CRR  from  Shear  Wave  Velocity 

Measurements  Along  with  Case  History  Data  Based  on  Revised  Values  of  MSF 
and  rj  Proposed  by  Idriss  (1999). 
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The  differences  between  the  MSF  and  relationships  proposed  by  Idriss  (1999)  and  the 
relationships  recommended  by  the  1996  NCEER  workshop  (Youd  et  al.,  1997)  are  significant 
at  magnitudes  less  than  about  7 (see  Fig.  2.4).  Figure  4.12  presents  two  liquefaction  resistance 
curves  for  earthquakes  with  magnitude  near  5.5  and  clean  soils  (FC  < 5 %).  The  upper  curve 
was  obtained  by  multiplying  values  of  CRR  defining  the  curve  for  FC  < 5%  in  Fig.  4.9  by  2.2, 
the  lower  MSF  recommended  by  the  1996  NCEER  workshop  for  magnitude  5.5  earthquakes 
(see  Eq.  2.13  with  n = -2.56).  The  lower  curve  was  obtained  by  multiplying  values  of  CRR 
defining  the  curve  for  FC  < 5 % in  Fig.  4.9  by  1.68,  the  MSF  proposed  by  Idriss  (1999)  for 
magnitude  5.5  earthquakes  (see  Eq.  2.15).  Also  shown  in  Fig.  4.12  are  the  available  case 
history  data  for  clean  sands  determined  using  average  stress  reduction  coefficients  proposed  by 
Seed  and  Idriss  (1971)  and  Idriss  (1998).  The  two  curves  in  Fig.  4.12  exhibit  differences  in 
CRR  of  about  0.02  at  V^j  = 100  m/s  and  0.1  at  = 200  m/s. 


4.3  RECOMMENDED  CRR-V^j  CURVES 

From  the  discussion  presented  in  Sections  4.1  and  4.2,  the  recommended  CRR-V^j 
curves  are  defined  by  Eqs.  2.13,  2.21,  and  4.1  with  a = 0.022,  b = 2.8,  and  n = -2.56.  The  value 
of  -2.56  for  n is  recommended  for  determining  magnitude  scaling  factors  because  it  provides 
more  conservative  CRR-V^j  curves  than  -3.3  for  magnitude  less  than  7.5.  While  the  magnitude 
scaling  factors  defined  by  Eq.  2.13  with  n = -2.56  provide  less  conservative  CRR-V^j  curves 
than  the  factors  proposed  by  Idriss  (1999)  for  magnitudes  less  than  7.5,  the  findings  of 
Ambrasey  (1988),  I.  M.  Idriss  (personal  communication  to  T.  L.  Youd,  1995),  Arango  (1996), 
Youd  and  Noble  (1997),  and  the  case  history  data  presented  in  this  report  (see  Figs.  4.2  and  4.3) 
support  their  use. 

The  recommended  CRR-V^j  curves  are  shown  in  Figs.  4.1  through  4.6  for  different 
magnitude  earthquakes.  The  curves  are  dashed  above  CRR  of  about  0.35  to  indicate  that  they 
are  based  on  limited  field  performance  data.  The  curves  do  not  extend  much  below  100  m/s, 
since  there  are  no  field  data  to  support  extending  them  to  the  origin. 
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Fig.  4.12  - Comparison  of  Liquefaction  Resistance  Curves  and  Case  History  Data 
for  Procedures  Recommended  by  the  1996  NCEER  Workshop 
(Youd  et  al.,  1997)  and  the  Revised  Procedures  Proposed  by  Idriss  (1999) 
for  Clean  Sands  and  Earthquakes  with  Magnitude  Near  5.5. 
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4.3.1  Correlations  Between  V^j  and  Penetration  Resistance 

By  plotting  values  with  equal  CRR,  on  can  obtain  Correlations  between  and 
corrected  penetration  resistance  from  the  recommended  CRR-V^j  relationships  given  in  Fig.  4.5 
and  the  1996  NCEER  workshop  recommended  SPT-  and  CPT-based  relationships  for 
magnitude  7.5  earthquakes  (Youd  et  al.,  1997). 

4.3.1. 1 Corrected  SPT  Blow  Count — Figure  4.13  presents  the  correlation  of  with 
(Nj)^q  for  clean  soils  (<  5 % fines)  based  on  the  recommended  CRR-V^j  and  CRR-(Nj)^q 
relationships.  Also  shown  are  the  field  data  and  mean  curve  for  sands  with  less  than  10  % non- 
plastic fines  from  Fig.  4.7.  The  correlation  derived  from  the  CRR  relationships  lies  between  the 
mean  and  the  mean  +15^^^  curves.  The  flatter  slope  below  (Nj)^g  of  6 exhibited  by  the  CRR- 
based  correlation  can  be  explained  by  different  assumed  minimal  values  of  CRR.  The  CRR-V^j 
relationship  for  magnitude  7.5  earthquakes  and  FC  <5  % shown  in  Fig.  4.5  provides  a CRR  of 
0.033  for  = 100  m/s,  the  lowest  value  shown  in  the  figure.  The  1996  NCEER  workshop 
recommended  a CRR  value  of  0.05  for  (Nj)^^  = 0.  The  difference  between  minimal  values  of 
CRR  is  small,  and  is  near  the  accuracy  of  both  procedures. 

The  CRR-based  curve  shown  in  Fig.  4.13,  along  with  the  plotted  field  data,  provide  a 
simple  method  of  comparing  the  and  (A^j^^-based  liquefaction  evaluation  procedures.  Both 
procedures  will  provide  similar  predictions  of  liquefaction  potential,  when  the  data  point  lies  on 
the  CRR-based  curve.  When  the  data  point  plots  below  the  CRR-based  curve,  the  -based 
liquefaction  evaluation  procedure  provides  the  more  conservative  prediction;  and  when  the  data 
point  plots  above  the  CRR-based  curve,  the  SPT-based  liquefaction  evaluation  procedure 
provides  the  more  conservative  prediction.  Since  most  of  the  data  points  shown  in  Fig.  4.13 
plot  below  the  CRR-based  curve,  the  V^^-based  procedure  provides  an  overall  more 
conservative  prediction  of  liquefaction  resistance  than  does  the  SPT-based  procedure  for  these 
sites. 


The  data  point  for  the  Treasure  Island  UM05  site,  which  incorrectly  lies  in  the  region  of 
no  liquefaction  shown  in  Fig.  4.9,  plots  just  below  the  CRR-based  curve,  as  shown  in  Fig.  4.13. 
Thus,  this  case  history  also  incorrectly  plots  in  the  region  of  no  liquefaction  on  the  SPT-based 
liquefaction  evaluation  chart.  Furthermore,  the  SPT-based  procedure  provides  a slightly  less 
conservative  prediction  of  liquefaction  resistance  than  the  shear-wave-based  procedure  for  this 
case  history. 

Although  the  CRR -based  curve  shown  in  Fig.  4.13  generally  trends  parallel  to  the  mean 
curve,  there  is  a small  hump  between  corrected  blow  counts  of  8 and  26.  This  hump  suggests 
that  either  the  CRR-V^j  relationship  is  more  conservative  or  the  CRR-(Nj)^q  relationship  is  less 
conservative  in  this  range. 
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4.3.1.2  Normalized  Cone  Tip  Resistance — Figure  4.14  presents  the  correlation  of  V^j 
with  for  clean  sands  with  median  grain  size,  between  0.25  mm  and  2.0  mm  based  on 
the  recommended  CRR-V^j  and  CRR-q^j^  relationships.  Also  shown  are  the  field  data  and 
mean  curve  for  clean  sands  with  less  than  10  % non-plastic  fines  from  Fig.  4.8.  The  correlation 
derived  from  the  CRR  relationships  lies  between  the  mean  and  the  mean  +15^^^  curves  for  > 
170  m/s,  indicating  that  the  V^^-based  procedure  provides  an  overall  more  conservative 
prediction  of  liquefaction  resistance  than  does  the  CPT-based  procedure  for  these  sites.  For 
< 170,  the  Ci?/?-based  correlation  lies  close  to  the  mean  curve,  indicating  that  both  procedure 
provide  an  overall  similar  prediction.  The  slope  of  the  CRR-ba.sed  correlation  below  q^j^  of  20 
may  be  explained  by  the  different  assumed  minimal  values  of  CRR. 

The  data  point  for  the  Treasure  Island  UM05  site,  which  incorrectly  lies  in  the  region  of 
no  liquefaction  shown  in  Fig.  4.9,  plots  on  the  CRR-bSiSed  curve,  as  shown  in  Fig.  4.14.  Thus, 
the  CPT-based  procedure  provides  a similar  prediction  of  no  liquefaction. 

4.3.2  Cementation,  Aging,  and  Above-the-Water-Table  Correction 


The  recommended  CRR-V^j  relationship  is  limited  to  the  characteristics  of  the  database, 
as  described  in  Chapter  3.  In  areas  of  cemented  and  aged  soils  (>  10  000  years),  and  soils 
above  the  ground  water  table  where  negative  pore-water  pressure  can  increase  the  value  of 
measured  in  seismic  tests,  a correction  factor  can  be  added  to  Eq.  2.21  as  follows: 


CRR  = 


a 


I 100  ) 


+b 


f 1 

^V*si-CVsi 


MSF 


(4.3) 


where 

C = a factor  to  correct  for  high  values  of  caused  by  cementation,  aging,  and 
negative  pore-water  pressures. 

Average  estimates  of  C range  from  0.6  to  0.8  based  on  the  penetration- correlations  for 
Pleistocene-age  soils  by  Rollins  et  al.  (1998a)  and  Ohta  and  Goto  (1978),  respectively. 
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Fig.  4.14  - Relationships  Between  and  V^j  for  Clean  Sands  Implied  by  the  Recommended 
CRR-V^j  Relationship  and  the  1996  NCEER  Workshop  Recommended  CRR-q^j^ 
Relationship  (Youd  et  al.,  1997)  with  Field  Data  for  Sands  with  Less  than  10  % 
Fines. 
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Figures  4.15  and  4.16  illustrate  a method  for  estimating  the  value  of  C using  SPT  and 
CPT  test  results,  respectively.  Shown  in  the  figures  are  the  V^^-penetration  correlation  curves 
for  clean  soils  and  silty  soils  implied  by  the  recommended  CRR-V^j  relationship  and  the  1996 
NCEER  workshop  recommended  CRR -penetration  relationships  (Youd  et  al.,  1997).  In  the 
example,  the  measured  values  of  V^j,  and  fines  content  are  220  m/s,  8,  55,  and 

10  %,  respectively.  The  -penetration  correlation  curves  suggest  a C value  of  0.74  for  these 
conditions. 

The  method  for  estimating  C described  above  assumes  that  the  strain  level  induced 
during  penetration  testing  is  the  same  strain  level  causing  liquefaction.  This  may  not  be  true 
since  pore-water  pressure  buildup  to  liquefaction  can  occur  at  medium  strains  in  several  loading 
cycles  (Dobry  et  al.,  1982;  Seed  et  al.,  1983). 

4.3.3  High  Overburden  Stress  Correction 

To  adjust  cyclic  resistance  ratios  where  overburden  stresses  are  much  greater  than 
100  kPa,  Seed  (1983)  developed  the  correction  factor  This  correction  is  applied  by: 

CRRj  = CRR  (4.4) 


where 


CRRj  = the  cyclic  resistance  ratio  at  the  reference  overburden  stress  (100  kPa). 

The  original  correction  factors  were  derived  from  isotropically  consolidated  cyclic 
triaxial  compression  test  results.  They  decreased  almost  linearly  with  effective  overburden 
pressure  from  a value  of  1.0  at  100  kPa  to  values  ranging  from  0.4  to  0.65  at  800  kPa. 
Modifications  to  these  factors  have  been  suggested  based  on  subsequent  analyses  of  additional 
cyclic  triaxial  and  constant-volume  cyclic  simple  shear  tests  (Vaid  et  al.,  1985;  Seed  and 
Harder,  1990;  Pillai  and  Byrne,  1994;  Vaid  and  Thomas,  1995;  Arango,  1996),  as  shown  in  Fig. 
4.17.  Figure  4.18  presents  the  relationship  of  minimal  values  of  recommended  by  the  1996 
NCEER  workshop  (Youd  et  al.,  1997;  Harder  and  Boulanger,  1997)  for  both  sands  and  gravels. 

Hynes  and  Olsen  (1998)  compiled  and  examined  a database  of  about  150  values  for 
gravels,  sands  and  silt-sand  mixtures  at  densities  ranging  from  very  loose  to  very  dense.  The 
database  consisted  of  both  undisturbed  and  reconstituted  specimens.  The  reconstituted 
specimens  were  pluviated  through  air,  pluviated  through  water,  and  constructed  in  layers  by 
moist  tamping  or  vibration.  Hynes  and  Olsen  compared  the  laboratory  estimates  of  CRR  at  one 
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Fig.  4.15  - Relationships  Between  (Nj)^q  and  V^j  Implied  by  the  Recommended  CRR-V^j 
Relationship  and  the  1996  NCEER  Workshop  Recommended  CRR-(Nj)^q 
Relationship  (Youd  et  al.,  1997)  with  an  Example  for  Determining  the 
Correction  Factor  C at  a Weakly  Cemented  Soil  Site. 
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Fig.  4.16  - Relationships  Between  and  Implied  by  the  Recommended  CRR-V^ 
Relationship  and  the  1996  NCEER  Workshop  Recommended  CRR-q^j^ 
Relationship  (Youd  et  al.,  1997)  with  an  Example  for  Determining  the 
Correction  Factor  C at  a Weakly  Cemented  Soil  Site. 
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Fig.  4.17  - Values  of  Determined  by  Various  Investigators  (after  Harder  and 
Boulanger,  1997). 


Fig.  4.18  - Minimum  Values  of  Recommended  by  the  1996  NCEER  Workshop  (after 
Youd  et  al.,  1997;  Harder  and  Boulanger,  1997). 
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atmosphere,  about  100  kPa,  with  CRRj  from  penetration  tests  in  the  field.  They  concluded  that 
method  of  deposition,  aging,  stress  history  and  density  strongly  influence  K^.  For  practical 
liquefaction  evaluation,  Hynes  and  Olsen  recommend  be  computed  as  (modified  for  a\  in 
kPa): 


(4.5) 


where 


/=0.8 

for  loose  soils. 

(4.6a) 

/=0.7 

for  moderately  dense  soils,  and 

(4.6b) 

/=0.6 

for  dense  or  slightly  overconsolidated  soils. 

(4.6c) 

For  very  dense  or  highly  overconsolidated  soils,  the  value  of  / may  be  less  than  0.6.  Figure 
4.19  presents  the  curves  defined  by  Eqs.  4.5  and  4.6  along  with  the  data  reported  by  Hynes  and 
Olsen.  The  recommendations  from  this  study  were  adopted  at  the  August  1998  MCEER 
workshop. 


Fig.  4.19  - Recommended  Values  of  Adopted  by  the  1998  MCEER  Workshop  (after  Hynes 
and  Olsen,  1998). 
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CHAPTER  5 


APPLICATION  OF  THE  PROCEDURE  FOR  EVALUATING 
LIQUEFACTION  RESISTANCE 


5.1  PROCEDURE  SUMMARY 

The  procedure  for  evaluating  liquefaction  resistance  outlined  in  the  previous  chapters  is 
summarized  in  the  following  steps: 

1.  From  available  subsurface  data,  develop  detailed  profiles  of  shear  wave 
velocity,  soil  type,  fines  content  (particles  less  than  75  jim),  and,  if  possible, 
soil  density  and  penetration  resistance. 

2.  Identify  the  depth  of  the  ground  water  table,  noting  any  seasonal  fluctuations 
and  artesian  pressures. 

3.  Calculate  the  total  and  effective  overburden  stresses  for  each  measurement 
depth  at  which  seismic  testing  has  been  performed. 

4.  Correct  the  shear  wave  velocity  measurements  to  the  reference  overburden 
stress  of  100  kPa,  as  shown  in  Eq.  3.2. 

5.  Determine  the  limiting  upper  value  of  overburden  stress-corrected  shear 
wave  velocity  for  liquefaction  occurrence,  V*s\ , for  each  measurement  depth 
using  Eq.  4.1.  If  the  fines  content  is  unknown,  assume  215  m/s  for  V51 . 

6.  Determine  the  value  of  the  correction  factor  C.  The  value  of  C can  be 
assumed  equal  to  1,  if  the  soil  to  be  evaluated  is  uncemented,  less  than 
10  000  years  old,  and  lies  below  the  ground  water  table  when  the  seismic 
tests  were  conducted.  If  the  soil  is  cemented,  more  than  10  000  years,  or  lies 
above  the  water  table,  the  value  of  C may  be  estimated  by  the  method 
described  in  Section  4.3.2.  If  the  soil  conditions  are  unknown  and 
penetration  data  are  not  available,  assume  0.6  for  C. 
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7.  Determine  the  design  earthquake  magnitude  and  the  expected  peak  horizontal 
ground  surface  acceleration. 

8.  For  each  measurement  depth  below  the  water  table,  calculate  the  cyclic  stress 
ratio,  as  shown  in  Eq.  3.1.  The  stress  reduction  coefficient  is  estimated  using 
Fig.  2.1. 

9.  Plot  values  of  and  CSR,  and  the  appropriate  liquefaction  resistance  curves 
using  Eqs.  2.13  and  4.3  with  a = 0.022,  b = 2.8,  and  n = -2.56.  If  the 
effective  overburden  stress  is  greater  than  100  kPa,  correct  curves  using  Eqs. 
4.5  and  4.6.  Liquefaction  is  predicted  at  the  site  if  the  data  points  plot  to  the 
left  of  the  respective  liquefaction  resistance  curve.  No  liquefaction  is 
predicted  if  the  data  points  plot  to  the  right  of  the  respective  curve. 


A common  way  to  quantify  the  potential  or  hazard  for  liquefaction  is  in  terms  of  a 
factor  of  safety.  The  factor  of  safety,  F5,  against  liquefaction  can  be  defined  by: 


FS  = 


CRR 

CSR 


(5.1) 


Liquefaction  is  predicted  to  occur  when  FS  is  less  than  1.  When  FS  is  greater  than  1, 
liquefaction  is  predicted  not  to  occur.  The  acceptable  value  of  FS  will  depend  on  several 
factors  including:  (1)  the  acceptable  level  of  risk  for  the  project,  (2)  the  extent  and  accuracy  of 
seismic  measurements,  (3)  the  availability  of  other  site  information,  and  (4)  the  conservatism  in 
determining  the  design  earthquake  magnitude  and  the  expected  peak  ground  acceleration. 

Thus,  the  following  step  may  be  added  to  the  nine-step  procedure  give  above— 

10.  Calculate  the  value  of  FS  for  each  measurement  depth  using  Eq.  5.1. 


5.2  CASE  STUDIES 

To  illustrate  the  application  of  the  procedure,  two  sites  shaken  by  the  1989  Loma  Prieta, 
California,  earthquake  (M^  = 7)  are  considered  below.  The  two  sites  are  Treasure  Island  Fire 
Station  and  Marina  District  School.  These  sites  are  discussed  earlier  to  illustrate  how  values  of 
and  CSR  are  determined  (see  Section  3.2). 


80 


5.2.1  Treasure  Island  Fire  Station 


Figures  5.1  and  5.2  present  the  liquefaction  evaluation  for  the  1989  earthquake  and  the 
crosshole  test  array  B1-B4  at  the  Treasure  Island  Fire  Station  site.  As  discussed  in  Section  3.2, 
values  of  and  CSR  are  calculated  assuming  densities  of  1.76  Mg/m^  above  the  water  table 
and  1.92  Mg/m^  below  the  water  table.  The  geometric  mean  of  the  two  peak  values  observed  in 
the  horizontal  ground  surface  acceleration  records  for  the  fire  station  and  the  1989  earthquake 
is  0.13  g.  Profiles  of  soil  type  and  fines  content  shown  in  Fig.  5.1  are  based  on  information 
provided  by  de  Alba  et  al.  (1994)  and  de  Alba  and  Faris  (1996).  The  value  of  C is  1,  since  the 
soil  to  be  evaluated  at  this  site  is  uncemented,  less  than  10  000  years  old,  and  lies  below  the 
ground  water  table. 


For  the  measurement  depth  of  4.6  m,  V^j  is  158  m/s  and  CSR  is  0.131  (see  calculations 
given  in  Eqs.  3.1  and  3.2).  The  average  fines  content  is  about  24  % (de  Alba  and  Faris,  1996). 
The  values  of  Vsi , CRR,  and  FS  are  calculated  by: 

Vli  = 215  - 0.5(FC-5)  = 215  - 0.5(24-5)  = 206  m/s  (5.2) 


and 


CRR  =U 


fCVsi 


r 


+b 


= |0.022 

= 0.119 


V 100 ; 
058^' 


1 


1 


Uoo 


{V*si-CVsi  V*sJ 

f 1 


4-2.8 


MSF 


1 


U06-158  206 


y Y2.56 

Ti, 


and 


FS 


CRR  _ 0.119 
CSR  ~ 0.131 


= 0.91 


(5.3) 


(5.4) 


Since  the  value  of  FS  is  less  than  1,  liquefaction  is  predicted. 
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Fig.  5.1  - Application  of  the  Recommended  Procedure  to  the  Treasure  Island  Fire  Sation  Site,  Crosshole  Test  Array  B1-B4  (Depth 
of  1.5  m to  14  m). 
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Fig.  5.2  - Recommended  Liquefaction  Assessment  Chart  for  Magnitude  7 Earthquakes 
with  Data  for  the  1989  Loma  Prieta  Earthquake  and  the  Treasure  Island  Fire 
Station  Site,  Crosshole  Test  Array  B1-B4  (Depths  of  1.5  m to  14  m).  Since  no 
Surface  Manifestations  of  Liquefaction  Occurred  at  the  Site  During  the 
Earthquake,  the  Data  Shown  are  for  a Non-liquefaction  Case  History. 
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Figure  5.1  shows  the  values  of  CSR,  CRR,  and  FS  for  crosshole  test  array  B1-B4  at 
the  measurement  depths  below  the  water  table.  Values  of  FS  are  less  than  1 for  the  depths  of  4 
m to  about  9 m.  Between  the  depths  of  4 m and  7 m,  the  sand  contains  non-plastic  fines  and  is 
considered  liquefiable  by  the  Chinese  criteria  (see  Section  1.1).  Between  the  depths  of  7 m and 
9 m,  the  soil  exhibits  plastic  characteristics  and  may  be  non-liquefiable  by  the  Chinese  criteria. 
Thus,  the  plotted  values  of  CRR  and  FS  for  this  clayey  sand  layer  are  shown  as  open  circles. 
The  layer  most  likely  to  liquefy,  or  the  critical  layer,  lies  between  the  depths  of  4 m and  7 m. 

Figures  5.3  and  5.4  present  the  liquefaction  evaluation  for  the  SASW  test  array. 
Locations  of  measurements  for  the  SASW  test  array  are  assumed  at  the  center  of  the  layer 
used  in  forward  modeling  of  surface  wave  measurements.  Values  of  FS  are  less  than  1 between 
the  depths  of  about  3.5  m and  1 1 m.  The  lowest  values  of  FS  in  the  non-plastic  soil  is  0.75  at  a 
depth  of  5.3  m.  This  FS  value  is  similar  to  the  lowest  FS  value  of  0.77  determined  from 
crosshole  measurements  in  the  critical  layer. 

Although  no  sand  boils  or  ground  cracks  occurred  at  the  fire  station  during  the  1989 
earthquake,  the  prediction  of  liquefaction  agrees  with  the  conclusion  stated  in  Section  3.2.1  that 
liquefaction  of  an  underlying  sand  cause  the  sudden  drop  in  the  acceleration  time  histories 
recorded  at  this  site  (de  Alba  et  al.,  1994).  It  is  possible  that  the  4 m thick  layer  capping  the 
site,  predicted  not  to  liquefy  (see  Figs.  5.1  and  5.3),  prevented  the  formation  of  sand  boils  at  the 
ground  surface  (Ishihara,  1985). 

5.2.2  Marina  District  School 

Figures  5.5  and  5.6  present  the  liquefaction  evaluation  for  the  Marina  District  School 
site  and  the  1989  Loma  Prieta  earthquake.  Values  of  V^j  and  CSR  are  calculated  assuming 
densities  of  1.76  Mg/m^  above  the  water  table  and  1.92  Mg/m^  below  the  water  table,  as 
described  in  Section  3.2.3.  They  are  plotted  in  Fig.  5.5  at  the  depths  midway  between  receiver 
locations.  The  profiles  of  soil  type  and  fines  content  shown  in  Fig.  5.5  are  based  on 
information  provided  by  Kayen  et  al.  (1990).  The  upper  7.6  m of  soil  at  the  site  consists  of 
sand  with  1 % to  8 % fines.  Since  the  sand  is  uncemented,  less  than  10  000  years  old,  and  lies 
below  the  ground  water  table,  the  value  of  C is  1. 

Calculated  values  of  FS  are  0.42,  0.90,  and  0.51  at  the  depths  of  3 m,  4 m,  and  6.7  m. 
The  clean  sand  is  underlain  by  a silty  clay  layer,  which  is  non-liquefiable  by  the  Chinese 
criteria  (see  Section  1.1).  Having  the  lowest  average  value  of  FS,  the  sand  fill  just  below  the 
water  table  between  the  depths  of  2.7  m and  4.4  m is  identified  as  the  critical  layer  that 
liquefied.  A prediction  of  liquefaction  agrees  with  the  observed  field  behavior. 
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Fig.  5.3  - Application  of  the  Recommended  Procedure  to  the  Treasure  Island  Fire  Sation  Site,  SASW  Test  Array  (Depths  of  2 m to 
13  m). 
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Fig.  5.4  - Recommended  Liquefaction  Assessment  Chart  for  Magnitude  7 Earthquakes 
with  Data  for  the  1989  Loma  Prieta  Earthquake  and  the  Treasure  Island  Fire 
Station  Site,  SASW  Test  Array  (Depths  of  2 m to  13  m).  Since  no  Surface 
Manifestations  of  Liquefaction  Occurred  at  the  Site  During  the  Earthquake, 
the  Data  Shown  are  for  a Non-liquefaction  Case  History. 
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Fig.  5.5  - Application  of  the  Recommended  Procedure  to  the  Marina  District  School  Site  (Depths  of  3 m to  10  m). 
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Fig.  5.6  - Recommended  Liquefaction  Assessment  Chart  for  Magnitude  7 Earthquakes 
with  Data  for  the  1989  Loma  Prieta  Earthquake  and  the  Marina  District  School 
Site  (Depths  of  3 m to  7 m).  Since  Surface  Manifestations  of  Liquefaction 
Occurred  Near  the  Site  During  the  Earthquake,  the  Data  Shown  are  for  a 
Liquefaction  Case  History. 
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CHAPTER  6 


SUMMARY  AND  RECOMMENDATIONS 


6.1  SUMMARY 

Presented  in  this  report  are  draft  guidelines  for  evaluating  liquefaction  resistance 
through  shear  wave  velocity,  V^,  measurements.  The  guidelines  were  written  in  cooperation 
with  industry,  researchers  and  practitioners,  and  evolved  from  workshops  held  in  1996  and 
1998.  They  include  the  development  of  a recommended  procedure  and  guidance  for  its  use. 

The  recommended  procedure  follows  the  general  format  of  the  simplified  penetration- 
based  procedure  originally  proposed  by  Seed  and  Idriss  (1971).  Cyclic  stress  ratios,  CSR,  are 
calculated  using  Eq.  2.1,  with  the  average  stress  reduction  coefficient  estimated  from  Fig.  2.1. 

measurements  are  corrected  for  overburden  stress  using  Eq.  2.7.  The  ten  main  steps  for 
applying  the  recommended  procedure  are  summarized  in  Section  5.1. 

Liquefaction  resistance  curves  are  based  on  a modified  relationship  between  overburden 
stress-corrected  shear  wave  velocity,  V^j,  and  cyclic  stress  ratio  for  constant  average  cyclic 
shear  strain  suggested  by  R.  Dobry  (personal  communication  to  R.  D.  Andrus,  1996).  The 
quadratic  relationship  proposed  by  Dobry  is  modified  so  that  it  is  asymptotic  to  some  limiting 
upper  value  of  V^j.  This  limit  is  related  to  the  tendency  of  dense  granular  soils  to  exhibit 
dilative  behavior  at  large  strains,  as  well  as  the  fact  that  dense  soils  expel  dramatically  less 
water  during  reconsolidation  than  loose  soils.  The  modified  relationship  is  given  by  Eq.  2.21. 
To  determine  the  unknown  parameters  of  Eq.  2.21,  liquefaction  and  non-liquefaction  case 
histories  from  26  earthquakes  and  more  than  70  measurement  sites  in  soils  ranging  from  clean 
fine  sand  to  sandy  gravel  with  cobbles  to  profiles  including  silty  clay  layers  are  analyzed. 
Penetration- V5  correlations  are  also  considered. 

Figures  4.1  through  4.6  present  the  recommended  evaluation  curves  for  uncemented, 
Holocene-age  soils  and  different  earthquake  magnitudes.  These  curves  are  defined  by  Eqs. 
2.13  and  2.21  with  a = 0.022,  b = 2.8,  Vsi  = 200  m/s  to  215  m/s  (depending  on  fines  content), 
and  n = -2.56.  These  values  were  selected  such  that  the  evaluation  curves  bounded  over  95  % 
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of  the  case  histories  where  liquefaction  occurred.  By  constructing  relationships  between 
and  penetration  resistance  from  the  recommended  evaluation  curves  and  plotting  available  in 
situ  test  data,  it  was  shown  that  the  -based  evaluation  curves  are  generally  more  conservative 
than  the  penetration-based  evaluation  curves.  Corrections  are  suggested  for  cemented  and  aged 
soils,  as  well  as  high  overburden  stress  conditions,  in  Sections  4.3.2  and  4.3.3. 

Caution  should  be  exercised  when  applying  the  procedure  to  sites  where  conditions  are 
different  from  the  database.  The  case  history  data  are  limited  to  relatively  level  ground  sites 
with  the  following  general  characteristics:  (1)  uncemented  soils  of  Holocene  age;  (2)  average 
depths  less  than  about  10  m;  and  (3)  ground  water  table  depths  between  0.5  m and  6 m.  All 
measurements  are  from  below  the  water  table.  About  three-quarters  of  the  case  history  data  are 
for  soils  with  fines  content  greater  than  5 %.  Almost  half  of  the  case  histories  are  for 
earthquakes  with  magnitudes  near  7. 

Three  concerns  when  using  shear  wave  velocity  as  an  indicator  of  liquefaction 
resistance  are  (1)  its  higher  sensitivity  (when  compared  with  the  penetration-based  methods)  to 
weak  interparticle  bonding,  (2)  the  lack  of  a physical  sample  for  identifying  non-liquefiable 
clayey  soils,  and  (3)  not  detecting  thin  liquefiable  strata  because  the  test  interval  is  too  large. 
The  preferred  practice  is  to  drill  sufficient  boreholes  and  conduct  sufficient  other  in  situ  tests  to 
detect  thin  liquefiable  strata,  to  identify  non-liquefiable  clay-rich  soils,  to  identify  silty  soils 
above  the  ground  water  table  that  might  have  lower  values  of  should  the  water  table  rise,  and 
to  detect  liquefiable  weakly  cemented  soils. 


6.2  FUTURE  STUDIES 

The  following  future  studies  are  recommended: 

1 . Additional  well-documented  case  histories  with  all  types  of  soil  that  have  and  have 
not  liquefied  during  earthquakes  should  be  compiled,  particularly  from  deeper  deposits  (depth  > 
8 m)  and  from  denser  soils  (U^  > 200  m/s)  shaken  by  stronger  ground  motions  > 0.4  g),  to 
further  validate  the  recommended  curves.  Also,  case  histories  from  lower  magnitude 
earthquakes  < 7)  may  improve  estimates  of  the  magnitude  scaling  factor. 

2.  Laboratory  and  field  studies  should  be  conducted  to  further  refine  estimates  of  V*si, 
the  limiting  value  of  V^j  for  liquefaction  occurrence.  For  example,  careful  laboratory  studies 
may  identify  more  clearly  the  influence  of  fines  content  and  particle  size  on  V*s\  • Additional 
careful  penetration- correlation  studies  may  also  help  refine  the  V*s\  estimates. 
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3.  Laboratory  studies  should  also  be  conducted  to  evaluate  the  implied  assumption 
observed  in  Fig.  4.9  that  at  low  values  of  (say  100  m/s)  liquefaction  resistance  is 
independent  of  fines  content. 

4.  Additional  work  is  needed  to  evaluate  the  significance  of  ignoring  soil  type  and 
horizontal  stress  in  the  overburden  correction. 

5.  Standard  test  procedures  exist  only  for  the  crosshole  test.  Standard  test  methods 
should  be  developed  for  the  other  in  situ  seismic  tests. 
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APPENDIX  A 


SYMBOLS  AND  NOTATION 


The  following  symbols  and  notation  are  used  in  this  report: 


A 

= 

parameter  that  depends  on  the  soil  structure; 

a 

= 

parameter  related  to  the  slope  of  the  CRR-V^j  curve; 

^max 

= 

peak  horizontal  ground  surface  acceleration; 

Bj.Bi 

= 

parameters  relating  and  penetration  resistance; 

b 

= 

parameter  related  to  the  slope  of  the  CRR-V^j  curve; 

C 

cementation,  aging,  and  negative  pore-water  pressures 
correction  factor; 

CPT 

= 

Cone  Penetration  Test; 

CRR 

= 

average  cyclic  resistance  ratio; 

CRR, 

= 

CRR  for  cyclic  triaxial  tests; 

CRRi 

= 

CRR  corrected  for  high  overburden  stress; 

C/^75 

= 

CRR  for  magnitude  7.5  earthquakes; 

CSR 

= 

cyclic  stress  ratio; 

DA 

= 

double-amplitude  axial  strain; 

Dr 

= 

relative  density; 

D50 

= 

median  grain  size  by  mass; 

exp 

= 

the  constant  e raised  to  the  power  of  a given  number; 

= 

age  and  soil  type  factors  for  correlating  and  Nf, 

FC 

= 

fines  content  (particles  smaller  than  75  |Lim); 

f 

= 

high  overburden  stress  exponent; 

= 

function  of  minimum  void  ratio; 

fiVaJ 

= 

function  of  average  peak  cyclic  shear  strain; 

G 

= 

shear  modulus; 

^max 

= 

small-strain  shear  modulus; 

Gn 

= 

^max  corrected  for  confining  stress  and  void  ratio; 

= 

secant  shear  modulus  at  7^^; 

g 

z= 

acceleration  of  gravity; 

K 

= 

coefficient  of  lateral  earth  pressure  at  rest; 

K 

= 

high  overburden  stress  correction  factor; 
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natural  logarithm  function; 
magnitude  scaling  factor; 
earthquake  moment  magnitude; 
stress  exponent; 

SPT  blow  count  in  Japanese  practice; 

SPT  energy-corrected  blow  count; 

SPT  energy-  and  overburden  stress-corrected  blow  count; 

magnitude  scaling  factor  exponent; 

reference  overburden  stress  (=  100  kPa); 

probability  of  liquefaction  occurrence 

factor  to  account  for  effects  of  multidirectional  shaking; 

shear  stress  reduction  coefficient; 

Spectral-Analysis-of-Surface-Waves; 

Seismic  Cone  Penetration  Test; 
sine  function; 

Standard  Penetration  Test; 
residual  standard  deviation; 
small-strain  shear  wave  velocity; 
overburden  stress-corrected  V^; 

limiting  upper  value  of  for  liquefaction  occurrence; 

mean  stress-corrected  V^; 

depth; 

function  of  depth; 
function  of  depth; 
average  peak  cyclic  shear  strain; 
mass  density  of  soil; 

cyclic  deviator  stress  in  cyclic  triaxial  tests; 
initial  effective  horizontal  confining  stress; 
mean  effective  confining  stress; 
initial  effective  confining  stress  in  cyclic  triaxial  tests; 
total  vertical  (or  overburden)  stress; 
initial  effective  vertical  (or  overburden)  stress; 
average  cyclic  equivalent  uniform  shear  stress  generated  by 
earthquake;  and 

maximum  cyclic  shear  stress  generated  by  earthquake. 
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APPENDIX  B 


GLOSSARY  OF  TERMS 

The  following  definitions  apply  to  this  report: 


Case  History 

An  earthquake  and  a test  array. 

Critical  Layer 

The  layer  of  non-plastic  soil  below  the  ground  water  table  where 
corrected  values  of  shear  wave  velocity  and  penetration  are  the 
least,  and  where  cyclic  stress  ratios  are  the  greatest. 

Liquefaction 

Occurrence 

Surface  manifestations  of  excess  pore-water  pressure  at  depth, 
such  as  sand  boils,  ground  cracks  and  fissures,  and  ground 
settlement. 

Moment  Magnitude 

An  earthquake  magnitude  scale  defined  in  terms  of  energy. 

Overburden  Stress- 

Corrected  Shear  Wave 
Velocity 

Shear  wave  velocity  measurement  corrected  to  a reference 
vertical  (or  overburden)  stress  of  100  kPa. 

Peak  Horizontal  Ground 

Surface  Acceleration 

The  peak  value  in  a horizontal  ground  surface  acceleration 
record  that  would  occur  at  the  site  in  the  absence  of  liquefaction 
or  excess  pore- water  pressures. 

Shear  Wave  Velocity 

The  velocity  of  a propagating  shear  wave  within  a material  with 
either  the  direction  of  wave  propagation  or  the  direction  of 
particle  motion  in  the  vertical  direction. 

Shear  Wave 

A body  wave  with  the  direction  of  particle  motion  transverse  to 
the  direction  of  wave  propagation. 
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Test  Array 


The  two  boreholes  used  for  crosshole  measurements,  the  borehole 
and  source  used  for  downhole  measurements,  the  cone  sounding 
and  source  used  for  seismic  cone  measurements,  the  borehole 
used  for  suspension  logging  measurements,  or  the  line  of 
receivers  used  for  Spectral- Analysis-of-Surf ace- Waves  (SASW) 
measurements. 
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APPENDIX  C 


SUMMARY  OF  CASE  fflSTORY  DATA 


Table  C.l  presents  a summary  of  case  history  data  described  in  Chapter  3,  and  used  in 
Chapter  4 to  establish  the  recommended  liquefaction  resistance  curves.  This  database  has  been 
expanded  and  modified  from  the  database  presented  by  Andrus  and  Stokoe  (1997).  Most  of  the 
modifications  are  minor  with  the  intent  to  have  the  data  conform  to  the  draft  guidelines 
presented  in  this  document.  The  major  modifications  are  based  on  new  information  or 
correction  of  an  error  in  calculations.  Some  case  histories  included  in  the  earlier  database  by 
Andrus  and  Stokoe  have  been  omitted  due  to  one  of  the  three  following  reasons:  (1)  The 
reported  average  downhole  measurement  is  for  a depth  interval  much  greater  than  the 
identified  critical  layer.  (2)  The  critical  layer  is  likely  older  than  10  000  years  and  contains 
carbonate.  (3)  The  location  of  the  critical  layer  or  field  behavior  is  uncertain.  References  for 
the  case  history  data  are  given  in  Table  3.1. 
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Table  C.l  - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.  1 - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  Vs  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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Table  C.l  - Summary  Information  for  -Based  Liquefaction  and  Non-Liquefaction  Case  Histories. 
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